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Dielectric Constant of Water from O° to 100° C 
C. G. Malmberg and A. A. Maryott 


An equal ratio arm, capacitance-conductance bridge, operated at frequencies below 100 
kilocycles per second, was used to measure the dielectric constant of water with an accuracy 
of better than 0.1 percent at 5-degree intervals over the range 0° to 100° C. At 25° C the 
dielectric constant was found to have the value 78.30, which is about 0.3 percent lower than 


that usually accepted. 
e= 87. 740 


with a maximum deviation of 0.01 


0.40008 + 9.398110) 2 


unit in dielectric constant. 


The data fit the equation 


—1.410(10-*) 8, 


The experimental method 


and sources of error are considered in some detail. 


1. Introduction 


Numerous determinations of the dielectric con- 
stant of water have been reported in the literature 
during the past half century. ! 
this work was not conducted with sufficient accuracy 
to provide data adequate for reference purposes. 
Several of the more recent investigators [1,2,3,4,5],? 
employing varied experimental techniques in which 
the accuracy was stated or has been inferred to be 
of the order of 0.1 percent, report values that are in 
close agreement at room temperature but not at 
higher and lower temperatures. Discrepancies 
amounting to a percent or more exist at 
temperatures. 

This investigation was undertaken to redetermine 
the dielectric constant of water over the range 0° to 
100° C with an accuracy better than 0.1 percent. A 
low-frequency bridge method with a Wagner earth 
to permit the use of three-terminal dielectric cells 
was employed. The various methods of measure- 
ment, generally some form of bridge or resonance 
method, are all more or less subject to errors assoc- 
iated with the residual impedances of the network 
that are enhanced when the medium possesses 
appreciable conductivity. The present assembly 
has the advantage of facilitating the minimization, 
ready control, and accurate evaluation of these 
residual errors, and of simplifying certain problems 
regarding the design of cells suitable for accurate 
absolute measurements. 

The values of dielectric constant reported in this 
investigation are significantly lower than previously 
assigned values [6], about 0.3 percent at 25°, but the 
variation of dielectric constant with temperature is 
substantially the same as reported by Wyman and 
Ingalls [5] except at the lowest temperatures. 


2. Apparatus 


The a-c bridge is of the equal-ratio arm, capaci- 
tance-conductance type, with a Wagner earth. 


The Wagner earth permits use of either two- or | 


three-terminal cells, and in conjunction with proper 


' For an extensive compilation and comparison of data prior to 1937, see N. E 
Yorsey, Properties of ordinary water substance, ACS Monograph, (Reinhold 
IF ublishing Co., New York, N. Y., 1940 

? Figures in brackets indicate the literature references at the end of this paper 
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However, much of | 


higher | 


shielding, either eliminates or fixes the influence of 
ground admittances and interbranch couplings. <A 
schematic diagram of the circuit is shown in figure 1, 
The basic network and its components, with the 
exception of the unknown arm and standard capaci- 
tor, are substantially the same as those described 
previously [7]. 

The unknown arm contains a shielded 3,000-ohm 
resistor, R, and a small capacitor, C, in pare llel with 
the guarded capacitance, C’,, of the test cell in order 
to provide convenient “zero” settings for the stand- 
ards, C,, and Ry, of the measuring arm. The capaci- 
tance of the test cell is measured in-terms of C, by 
connecting the cell in parellel first with C and then 
with 7, of the Wagner earth by means of the switch, 
K. Jn the latter case the cell is effectively removed 
from the main bridge circuit because its earth- 
potential electrode is connected to ground, and 
admittances from any corner of the bridge to ground 
do not affect the conditions of balance. 





Figure 1. Schematic diagram of the bridge circuit. 











The three cells used were designed for “absolute” 
measurements, the insulating supports being located 
where they should not affeet the desired interelec- 
trode capacitances. Intercomparisons between the 
three cells, using media of various dielectric con- 
stants, indicated that these capacitances varied 
linearly with the dielectric constant of the medium. 
Because accurate measurements become increas- 
ingly difficult as the conductivity of the medium 
increases, the materials of construction of the two 
cells, A and #2, used with water, were selected to 
minimize electrolytic contamination. Cell C was 
used only for intereomparisons with liquids of lower 
dielectric constant and conductivity. 

Cell A is a fixed, three-terminal, evlindrical capaci- 
tor. A schematic drawing is shown in figure 2. 
The material is copper, with the electrode 
surfaces coated with pure tin. The centre, or high- 
potential, electrede, /7, is supported by the grounded 
guard electrode, G. Electrical insulation and a 
liquid-tight seal are provided by 0.01-in. Teflon 
gaskets on each side of the silvered copper spacer, 
F, and by a Bakelite bushing on the upper side. 
The high-potential terminal is electrostatically iso- 
lated from the case by the removable shield, S. 
The guard electrode is supported on the outer case, 
which, in part, serves as the earth-potential electrode 
D, by means of a Pyrex-glass ring and 0.01-inch 
Teflon gaskets, and held by six bolts insulated from 
the guard electrode by Bakelite bushings. To pre- 
vent the trapping of gas in the annular space between 
the guard and the earth potential electrodes when 
the cell is being filled with liquid, four 1-mm holes 
were drilled through G at level &. A Pyrex-glass 
reservoir is attached to a_ silvered copper tube 
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Figure 2 Sectional view of cell A. 
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extending from the bottom of the cell, as shown ay 
the left of figure 3, to facilitate the introduction: and 
manipulation of fluids. Coaxial cable, with the 
outer conductor grounded, connects the cell to the 
appropriate bridge terminals. Coupling 
the exposed, earth-potential electrode and othe: com- 
ponents of the network is prevented by proper 
shielding. As lead, shield, and guard capacitances 
are all effective directly to ground, only the capaci- 
tance between the test electrodes, D)) and IT, is 
measured. This capacitance in air is about 34 yuf 

Cell #, shown at the right in figure 3, 
differential type. It was constructed from a Gen. 
eral Radio type 568 D variable air capacitor modified 
in several respects. The number of plates was re. 
duced to give a replaceable capacitance of about 28 
uul, and all metal surfaces were tinned. The steatite 
insulators were replaced with Pyrex glass and_par- 
tially shielded to minimize any effect of the insulation 
on the differential capacitance. The unit was 
mounted in a gold-plated brass case. The rotor 
control handle is variable between fixed — stops 
mounted on the cover. The bigh-potential stator 
plates bracket the earth-potential rotor plates when 
interleaved to avoid varving the stray fields through 
the insulating supports. Both are isolated 
the case except when the cell is two- 
terminal capacitor. Then the rotor plates and the 
case are connected together. A comparison (table 1) 
of the values of dielectric constant obtained for a 


used as a 


given medium when the cell is used as both a two- 
and a three-terminal capacitor provides a check of 
the linearity of its replaceable interval of capacitance 
with respect to dieleetrie constant. Any nonlinear 


effeets associated with the insulating supports should 
differ in the two cases. 





Figure 3. Cell A (left) and cell B (right) 
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Cell - also of the differential type, and similar 
in des io cell B. It consists of a General Radio 
type Joo) variable air capacitor, mounted in a 
brass tainer, and has a replaceable capacitance 
of ab 500 wef. Both electrodes are insulated 
from case to permit use as either a three- or a 
two-tc! ‘inal capacitor. 

The standard capacitor, C,, is a General Radio 
type 722-D variable precision capacitor with the 
follow: modifications: The original wood ease, 
lined l) sheet copper, was replaced by a metal 
ease to improve the stability of capacitance; in 


addition, an external worm and gear drive was 
mounted on the panel to expand the seale and to 


allow a precision of setting of 0.01 uuf on the 100 to 
1.100 wuf range. The capacity of this interval was 


assumed to be exact and was calibrated internally 
in the conventional stepwise manner, using the re- 
placeable air capacitances of cells A and B and other 
larger and smaller fixed increments. Calibrations 
made several times during the course of the measure- 
ments on water and with the different increments 
were consistent to 0.05 yuf or better throughout the 
range of the precision capacitor. The total backlash 
associated with the setting of this capacitor amounts 
to about 0.005 wuf on the 100 to 1,100 pul range and 
about one-tenth of this value on the 25 to 110 wuf 
range. Mica capacitors, General Radio type 505, 
were plugged to the terminals of the precision 
capacitor to extend its range as needed. These were 
calibrated in terms of the standard capacitor at the 
required frequencies, and over a period of a year 
their capacitances were found to remain constant to 
better than 0.01 pereent. 

The apparatus was contained in a constant-tem- 
perature room maintained at 23°C. The cells were 
placed in a grounded oil bath controlled to better than 
0.01 deg and the temperatures determined with a 
platinum resistance thermometer. Cell A was im- 
mersed to within L inch of the top and the exposed 
portion provided with adequate thermal shielding. 
Experimental evidence showed that without such 
shielding, significant thermal gradients existed within 
the cell at the higher temperatures. Cells B and ¢ 
were immersed so that the oil just covered the lids 
and were used only for measurements at 25° C. 


3. Water 

Distilled water was obtained as required from the 
regular supply piped to the laboratory. This supply 
is free of ammonia and has an average conductivity 
of 0.5 micromho/em at 25°, due mainly to dissolved 
COQ,. Each sample was boiled in the test cell prior 
to measurement. This procedure was effective in 
removing dissolved air in order both to eliminate the 
tendency for bubbles to form between the electrodes 
and to reduce the COQ, content, thereby decreasing 
An atmosphere of hydrogen was 
thereafter maintained in the cell. After thorough 
cleaning and leaching of the cell, it was possible to 
keep the conductivity of the water in the range 0.1 


) to 1.7 micromho/em at all temperatures and to avoid 


rapid drifts, which would interfere with the precision 
of capacitance balance. 


4. Procedure and Errors 


If (’. is the true interelectrode capacitance of the 
fixed, guarded cell, A, when filled with a medium ot 
dielectric constant, e,, then 


fm (1) 


where (, is the capacitance in vacuum, C, the capaci- 
tance in air, and e, the dielectric constant of air. If 
the residual impedances associated with the bridge 
network are negligibly small, the interelectrode 
capacitance is obtained directly as the difference in 
capacitance of the standard capacitor, (, jal 
where ©, and ©, are the capacitances measured 
with the cell connected to the appropriate bridge 
terminals and with the earth-potential electrode 
connected directly to ground, respectively. In prac- 
tice, errors that may seriously limit the accuracy of 
the measurements of dielectric constant can arise 
from the effects of residual impedances associated 
with the leads and components of the bridge circuit. 
Although these errors can be minimized by careful 
design of the network and the selection of com- 
ponents, by optimum choice of frequency and cell 
capacitance, and by avoiding excessively high con- 
ductivities resulting from electrolytic contamination 
of the medium, it is often not feasible to reduce these 
errors to negligible proportions, and procedures 
must be devised for their evaluation. 

If the admittances between each test electrode 
and the guard electrode of the cell are neglected for 
the moment, the actual network can be simulated 
by that shown in figure 4. The effective capacitances 
of the ratio arms are represented by C, and (), the 
self-inductance of the standard condenser and its 
leads by a the effective series inductance of the 
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Figure 4. Schematic circuit showing the distribution of 
residual impedances in the bridge network. 








resistance box by J,, the cell resistance by R,, and 
the inductance and resistance of the leads to the 
cell by J and r. The residual impedances associated 
with the fixed 3,000-ohm resistor, R, and smali con- 
denser, C, of the unknown arm are neglected, as they 
do not alter the error analysis to follow. The cell 
capacitance is given by 


where, neglecting terms of higher order, 


Ce (C} C)(R, R, R, R}). 
Cyl) R2+-1/ Re —L,/R3, 
OF 2rC,/R,, 
Ca= ln (C2— CL) — wl C2? So (In —D) C2. 
C,, the capacitance due electrode polarization, 


varies approximately as w~", where @ is the angular 
frequency and n» generally lies between 1.5 and 2 but 
is virtually constant for a given experiment and 
temperature. The unprimed symbols refer to the 
case where the cell is connected in the normal man- 
ner, and the primed symbols to the case where the 
earth-potential electrode is connected directly to 
ground. 

The errors attributable to Cy, and C,, being fune- 
tions of the frequeney, were eliminated in the fol- 
lowing manner. In each instance values of C,, were 
obtained at 3, 6, 12, 24, 48, and 96 ke. A representa- 
tive plot of such data is shown by the solid line of 
figure 5. At the lower frequencies the variation of 
C', is due predominately to electrode polarization, 
and the error associated with the lead inductances 
is negligible or greatly minimized, whereas the con- 
Verse is true at “the higher yt g° neies. The unbal- 
ance in lead inductances, /,,—/, which amounted to 
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FREQUENCY, kc 


Capacitance Cy» plotted on an inverse frequency 
scale. 


Ficure 5. 


m after correction for Cg. 


_ measured values 


| 


was determined experimentally 
appropriate correction made, 


1.2 wh, and the 
as indicated by the 
broken line of figure 5. Extrapolation of this latter 
curve to infinite frequene y then gives the cell capaei- 
tance, corrected for the frequency-dependent errors. 
Determinations of /,,—/ were made both during the 
normal course of the measurements on water at the 
lowest temperatures, where electrode polarization was 
negligible at all but the lowest frequencies, and by 
plug rging sufficient capacitance in the form of the 
mica standards across the empty cell and standard 
condenser 
capacitance balance with frequency. The uncer. 
tainty in the value of the dielectric constant of water 
associated with these corrections is definitely 
than 0.01 unit. 

Of the remaining terms in eq (2), C, is not signif. 
icant in the present case as long as good electric 
contacts are maintained. The sum of the errors 
(',+C>, which depend upon the resistance balance, 
was obtained by a modification of the procedure 
previously described [7]. The dielectric cell was 
replaced by a conductivity cell of negligible capaci- 
tance, with J kept very nearly the same. The 
resistance of the cell filled with distilled water was 
then decreased in small steps by the addition of a 
solution of KCl. The change in C,, required to 
rebalance the bridge was observed at about 10-ohm 
intervals. The entire range of Ry required for the 
measurements on water, namely, 1,000 to 3,000 
ohms, was calibrated in this manner, using a fre- 
quency sufficiently high to avoid any effects of 
electrode polarization. The range 1,000 to 2,600 
ohms was also calibrated directly, using cell A. 
In this case the resistance was changed by the addi- 


less 


tion or removal of small amounts of CO... Both 
calibrations were in satisfactory agreement. The 
correction for (,+(C,, which, in the extreme case, 


amounted to about 0.07 unit in dielectric constant, 
should not lead to an uncertainty greater than 
0.005 unit. 

Determination of the replaceable capacitance of 
the cell when filled with dry air, C,, was made as 
part of the procedure employed for calibrating the 
standard capacitor. Thus, C, was equal to the total 
interval of capacitance that could be covered in 
integral steps of C, divided by the number of steps. 
Four determinations of C,, made at various times 
during the measurements on water, gave an average 
deviation of 0.007 percent. The variation of C, 
with temperature was obtained from measurements 
at a number of temperatures in the range 0 to 100°, 
using the low range of the standard condenser, where 
a precision of setting of 0.001 uuf was possible. 
When C, is corrected to vacuum, using adequate 
values for the dielectric constant of air [8], the 
vacuum capacitance, C,, is given by 

(C,) -=33.76,[1 + 1.8)(10~°) (t— 25) Juul. 

The change with temperature is -_ nearly the 
same as predicted from the linear coefficient of 
expansion of copper. 


to cause a readily measurable change of! 
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ninations ¥ the dielectric constant of water 


Jeter 
eo. made in cell A at 5-deg intervals, using two 
or mor sample s of aaa, at each — As 
a further check on the correction for C,+C», the use 
of samples of corresponding conductiv ity was avoid- 
ed. Without exception, the values of dielectric 
constant obtained at each temperature agreed to 


0.01 unit or better with respect to their mean. 

As an over-all check on the accuracy of the meas- 
urements with cell A, some additional experiments 
were made at 25° C, using auxiliary cells B and C 
both as two- and as three-terminal cells. The na- 
ture of the residual errors and the procedures used 
in their evaluation are similar to those discussed 
above for cell . However, as these cells are of the 
differential type, they require twice as many meas- 
urements of capacitance and associated corrections. 
The dielectric constant of water was redetermined 
in cell B, whereas ethylene chloride and 2-ethyl- 
hexanol were measured in cell C. The latter ma- 
terials were of technical grade. The degree of purity 





was immaterial as ‘a ‘sample was used for inter- | 


comparison in cell A. The values of dielectric con- 
stant obtained with the three cells are listed in table 
1. The differences are quite negligible except in 
the case of water, where cell y gives a value about 
0.03 percent larger than cell . However, as the 
experimental uncertainties .* hosed with cell 5 are 
at least twice those of cell A, this difference is hardly 
significant. 


TABLE 1 Intercomparison of values of dielectric constants 
measured at 25° C with the diffe rent cells 
Dielectric constant 
Cell 

Water CcHysO* C,HyCh* 

A 78. 30 6. 831 10. 376 
,f2-terminal 78. 32 
"\3-terminal 78. 32; 

f2-terminal 6. 831, 10. 377 

(3-terminal 6.831, 10. 378 


® Technical grade, unpurified 


It was assumed in the previous analysis that the 
test cell and its leads are adequately represented by 
the series-parallel combination of r, 1, and C,, and 
R., shown in figure 4, for which the equivalent 


parallel capacitance is given to a sufficient approxi- | 


mation by 


+ wIC',—2r/R,)—lU/R?. 3) 
A more rigorous consideration of its properties as a 
three-terminal impedance with associated leads re- 
sults in the equivalent circuit of figure 6. D, B, and 
@ are the points of connection to the bridge and to 
ground, and D’, B’, and G@’ are the corresponding 
electrodes of the cell. Analysis of this circuit was 
made by using the star -mesh transformation [9]. 
The equivalent parallel capacitance is given to a 


sufficient approximation by 
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Figure 6. Schematic circuit for the three-terminal cell. 
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It is evident from this relation that significant lead 
impedances can seriously limit the utility of a three- 
terminal cell unless the admittances from each of 
the test electrodes to the guard electrode are kept 
small compared to the direct admittance between 
the test electrodes. The procedures described pre- 
viously for evaluating the residual errors are not en- 
tirely valid when the effective capacitance is caleu- 
lated by eq (4) rather than eq (3). In the present 
measurements, however, where the admittances to 
the guard electrode were comparatively small, com- 
putations indicate that the use of eq (3) should not 
lead to an error in dielectric constant greater than 
0.01 unit in any case. In fact, if any difficulty had 
arisen from this simplification, it would have be- 
come evident at various stages during the course of 
the measurements, as, for example, in a discrepancy 
in the calibrations for C,+C, when the conductivity 
and dielectric cells were used. 


5. Results 


The values of the dielectric constant of water 
measured over the range 0.1° to 99° C fit the equation 


(5) 
The 


«= 87.740 —0.4008t+ 9.398(107-*)#—1.410(10~*)8, 


where ¢ is the temperature in degrees Celsius. 








values computed from this equation are listed in 
table 2, together with the differences between these 
values and those actually measured at each of the 


specified temperatures. The maximum difference 
does not exceed 0.01 unit. The last two columns 
list the values of de/dt and the temperature coefficient 
(l/e)(de dt). It is interesting to note that the tem- 
perature coefficient ts remarkably constant over 
the entire range. In fact, the alternative relation 

logye= 1.94315 —0.00197208, (6) 
fits the data nearly as well as eq (5). In this case 
the maximum deviation 0.02 and the average 
deviation 0.01 unit in dielectric constant. 

It is estimated that the over-all accuracy of the 
tabulated values of dielectric constant is +0.05 unit 
or better, and that the temperature coefficient. is 
determined to 1 percent or better at all temperatures. 


Is 


TABLE 2 Dielectric constant of water and related data at 
various temperatures 
, . 1 de 
eed €eq e obs ded + 
c 
uv S7. 740) . 4001 1 ox 10 
it S7. 700 +0) OM 
’ S5. 7H3 ol S008 4. SAT 
iL S34. S32 wre ss17 } $ 
| 1. O46 4 STU i ‘ 
1) sO. 103 OOS so 4.547 
~ 78. 304 0s $557 4.545 
ah 7H. 46 O08 BATS $. 534 
$5 74.828 (we i3U 4. 537 
Ww 73. 151 ore $317 4. 4 
} 71.512 (4H $241 4. 532 
4) wy uo Lo 4107 1. Sau 
i Os. 34 mm sas +. 520 
mo) wi Sl ia yr 1 » 
ih 65. 319 oo Us 4. 52s 
70 3. SS7 mr 2s $520 
rt 02.427 ~ 1 Pan al) + 531 
st) OL. O27 00s J7Hs 4. SBS 
s 54. 159 w] 2704 1. 41 
") “S31 OO M2 4 547 
‘ 7. O07 | 2597 $. 55 
" wT IS 
loo 7 2544 4. tit 


6. Discussion 


A comparison of the present data with those re- 
ported, over a wide range of temperature, by certain 
investigators in the last quarter century is shown in 
figure 7. The differences between the values of di- 
electric constant obtained in the earlier work and the 
present investigation are plotted as a function of the 
temperature. Smooth curves are shown where the 
data were fitted to empirical relations. With re- 
spect to the variation of dielectric constant with 
temperature, only the data of Wyman and Ingalls [5], 
which were obtained in terms of Wyman’s [1] value 
at 25° C, are in satisfactory agreement with the 
present measurements over an extended interval of 
temperature. Although their values are consistently 


higher, the temperature coefficient derived from the 
two sets of data 
all temperatures 


agrees to better than 1 percent at 
above 20° ©, 
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Figure 7. Comparison of the differences in the values of the 


dielectric constant of water obtained by various workers 


Ae =e literature—e present work 


Drake, Pierce, and Dow (1930 
Wyman and Ingalls (193s 
Wyman (1930 

Akerlof and Oshry (1950 11 
Lattey, Gatty, and Davies (1931 


Albright (1937 


io) Albright and Gosting (1946 
lr. T. Jones and R. M. Davies (19389 
Mean of 17 values selected from literature, Lattey, Gatty, and 
Davies (1931 


Various reported values of the dielectric constant 
at 25° C are compared in table 3. The value of 
Lattev, Gatty, and Davies [10] resulted from a sur- 
vey of the literature prior to 1930. Although this 


value is closest to that obtained here, the agreement) 
is not significant because of the wide seatter in the} 


individual results (mean deviation 0.4 unit). The 
remaining values, each obtained by a different. ex- 
perimental procedure and considered to have an 
accuracy of the order of 0.1 percent, are roughly 0.3 
percent higher than that reported in this investiga- 
tion. The agreement between these earlier 
values at 25° C would appear to be rather fortuitious 
in view of the sizable discrepancies existing at other 
temperatures, as shown by figure 7. 


close 


TABLE 3. Comparison of values of the dielectric constant of 
water at 25° ¢ 

e NIethod Frequet Auth 

78. 25 Mean of 17 reported Latteyv, Gatt Da 
values prior to 1430 vie 1431 

‘Ts. 54 Resonators suspended fto sl Nike Wyman (1930 
in mediun 

78. 57 Standing waves in co l2Zto 77 Me Drake, Pierce, Dow 
iXical Linn 1u30) 

TS 48 Bridges ATO ke Albright (1937 

‘TS. 40 Voltage n OTLATECE O70 ke Tones and Davies 

1u3y 


7S. 30 Bridge 3 to OW ke Present work 


* Values apparently relative to air rather than to vacuum 
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Alt! the work of Jones and Davies [4] was 
of hig! cision and conducted with unusual care 
and xt on to possible sources of error, the method 
was a see ative one in which benzene was em- 
ploved ‘primary standard. Consequently, anv 
error 1 cominan value for the dielectric constant 
of benz would result in a percentage error nearly 
twice as large in the value for water. Furthermore, 
direet comparisons of water and benzene in the same 
eell were vot feasible, so that additional cells and 


liquids o! intermediate dielectric constant were re- 
quired. ‘The number of operations upon which the 
value for water was dependent was thereby increased 
and the over-all uncertainties of the me thod somewhat 
enhanced. In view of these limitations the difference 
between this and the present work does not appear 
to be excessively large. In fact, 
and Davies are corrected for the 
recently recommended [6], this 
duced to about 0.1 percent. 

The remaining values in table 3 were obtained by 
methods essentially absolute in that no standard 
medium other than air was needed. Because of the 
lack of experimental detail pertinent to the possible 
sources of error in the report of Albright [3], no de- 
tailed comment on this work is possible. The ex- 
perimental methods of Wyman [1] and of Drake, 
Pierce, and Dow [2] were comparatively direct and 
simple in principle. Although the experimental pre- 
cision was less than that obtainable with conven- 
tional bridge and resonance methods, this disad- 
vantage was seemingly offset by the reduction or 
avoidance of troublesome residual errors. In both 
eases the value reported at 25° C was the mean of a 
number of values obtained at various frequencies 
which had a spread of more than 0.3 percent. 

The method of Wyman involves the determination 
of the resonant frequencies of fixed, metallic 
nators suspended in the medium. Seven resonators 
of several different designs and of varying frequencies 
were emploved. Exclusive of the value obtained 
with the resonator of lowest frequency, which was 
about 5 percent too high, the values of dielectric 
constant, when corrected to 25° C, ranged from 78.42 
to 78.70. It was assumed that these resonators 
behaved as idealized lumped circuits of inductance 
Land capacitance (, so that the resonant frequeney 


value of benzene 
difference is re- 


reso- 


f=1/[20(LO)'7].) Then «=/f, where fy is the 
resonant frequency in vacuum, and f is the corre- 
sponding frequency in the medium. The. strict 


validity of this relation is questionable because of the 
effect of conductance of the medium on the properties 
of resonant circuits and the dependence of the in- 
ductance and resistance of metallic conductors upon 
frequency through the “skin effeet.””. A change in 
inductance due to the skin effeet would lead to a high 
value for the dielectric constant. Although the mag- 
nitude of this error is difficult to evaluate without 
more specific details regarding the construction of the 


if the data of Jones | 


resonators, estimates based upon simple U-loops, 
which simulate the form of inductance in some of thu 
resonators, indicate that an error of the order of 
several tenths of 1 percent is plausible. From a 
study of the behavior of several resonators when the 
ionic conductivity of the water was increased by the 
addition of small amounts of KCl, Wyman showed 
that the resulting error also gave high values for 
the dielectric constant, but concluded that the error 
was not significant in the case of pure water. How- 
ever, estimates based upon the data presented, after 
making allowance for the small but significant con- 
tribution to the conductivity arising from dielectric 
indicate a remaining error of the order of 0.1 
percent in most instances. 

The procedure of Drake, Pierce, and Dow in- 
volved, in essence, the determination of the half 
wavelength of the standing waves set up by a source 
of known frequency in a 4-m, vertical, coaxial pipe 


loss, 


filled with water. Except for factors computed 
to be of negligible sper the dielectric con- 
stant was ae by «=c’ (f?Xm”) , Where f is the 


frequency, A» the wavele ngth in the medium, and ¢ 
the velocity of light. The precision obtained in 
determining the values of dielectric constant was 
governed primarily by the limitations in measuring 
half wavelength and temperature, about 0.2 mm 
and 0.2° C, respectively, and varied from 0.3 to 0.15 
percent, depending upon the frequency. Inasmuch 
as the presence of systematic errors of comparable 
magnitude would hardly be detected, the difference 
noted between their value and the present value is 
not unreasonable. 

The magnitude of the disagreement shown in figure 
7 indicates that, in general, the accuracy of these 
data can hardly approach 0.1 percent. This is not 
too surprising as the reproducibility obtained was 
only rarely better than 0.1 percent. 

For the present investigation, the uncertainty due 
to any known source of error seems to be limited to 
-0.01 unit. Summation, without regard to sign, 
of the individual uncertainties associated with the 
recognized sources of error, leads to a possible un- 
certainty of +0.05 in the values of € or +0.1 percent 
at the highest temperature. The over-all repro- 
ducibility of separate determinations at any single 
temperature, ranging in number from 2 to 10, was 
such that the largest deviation from the mean for 
any given temperature was never worse than +0.015 
unit, which occurred with both cells at 25° C 
(4+6=10 observations); and the average of all 69 
deviations from the individual means for 21 temper- 
atures was +0.005 unit. This agreement indicates 
the of significant cumulative error and 
should provide a reasonable basis for an estimation 
the accuracy attained. As a consequence, the 
assignment of an accuracy of +0.05 unit to all these 
data seems entirely justified. 


absence 


of 
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A Calorimeter for Measuring the Power in 
a High-Energy X-ray Beam* 


John McElhinney,' Bernard Zendle, and Steve Domen 


The design and calibration of a calorimeter to measure the power in X-ray beams 


having peak energies between 1 and 180 million electron volts are described. 


The calorim- 


eter included two thermally balanced lead cylinders, 4 centimeters in diameter by 7.5 


centimeters long, one irradiated by an X-ray beam, 


absorb almost completely the X-ray beam. 


embedded thermistors. 


The lead cylinder was large enough to 


The absorbed energy resulted in an unbalance 
of temperature of the two cylinders, which was measured by the change in resistance of 


Calibration of the calorimeter consisted in observing the temperature rise due to a measured 


quantity of electric energy dissipated in the same cylinder. 


The results are given for five 


calibration runs, each using about 70 microwatts of power for approximately 20 minutes. 


The probable error of the mean was about +1 percent. 


Separate reports of measurements 


of X-ray-beam powers at 1.4 and 36 million electron volts are in preparation. 


1. Introduction 


The interpretation of physical and medical experi- 
ments with high-energy X-rays has been greatly 
hampered by the uncertainties in the measurement of 
X-ray-beam power. The conventional method of 
determining beam power has been to measure the 
ionization current in the gas of an ionization chamber 
placed in the X-ray beam. To calculate the relation- 
ship between the beam power and the ionization 
current, one must make various assumptions con- 
cerning the spectrum of the incident photons, the 
production of electrons and positrons by the photons 
in the chamber walls, the values of the stopping power 
for these electrons in the wall and the gas, and the 
energy required for an electron or positron to produce 
an ion pair in the gas of the ionization chamber. 
Because these assumptions and quantities are diffi- 
cult to evaluate, the calculated response of an 
ionization chamber to X-rays may contain unknown 
errors of the order of 5 to 10 percent. 

In order to avoid these uncertainties, apparatus 
was constructed at the National Bureau of Standards 
for measuring the X-ray-beam power calorimetri- 
cally. Such measurements were made by comparing 
the temperature rise of a lead cylinder irradiated by 
the X-ray beam with the temperature rise of the 
same cylinder heated by a measured quantity of 
electric energy. After making suitable corrections 
for the small proportion of X-ray energy that was 
not absorbed, an electrical equivalent of X-ray 
energy could be obtained. 

The first published measurements of this type 
were made by Laughlin et al. [1]. His later publi- 
cation [2] described an improved calorimeter used to 
measure beam power from a 22.5-Mev betatron and a 
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400-kev * X-ray unit. Laughlin’s calorimeter con- 
sisted of twin lead cylinders, one of which was 
irradiated. The resulting temperature difference 
between the two cylinders was used to determine the 
total energy absorbed. 

Calorimetric measurements of X-ray-beam power 
were also made by Edwards and Kerst [3] for X-ray 
peak energies ranging between 150 and 300 Mev. 
Their calorimeter was similar to Laughlin’s in 
design, but, due to the higher energy X-rays, more 
consideration had to be given to the power escaping 
from the cylinders. 

One of the projects of the Betatron Section at the 
National Bureau of Standards was to develop an 
instrument to measure the power in X-ray beams 
with peak energies between 1 and 180 Mev. This 
report describes the construction and operation of a 
calorimeter designed for this purpose. X-ray-beam 
power measurements at 1.4 and 36 Mev with this 
calorimeter have been made, and the results are now 
being prepared. 


2. General Method of Operation of the 


Calorimeter 


The calorimeter consisted, essentially, of two 
similar lead cylinders, thermally insulated from their 
surroundings and from each other. The X-ray beam, 
emerging from an aperture in the shielding wall 
around the X-ray source, was directed axially into 
the end of one of the lead cylinders. The X-ray 
energy absorbed in the cylinder was converted into 
heat energy, resulting in a rise in temperature of the 
lead cylinder. The temperature-sensing elements 
were thermistors embedded in the two cylinders and 
connected in a Wheatstone-bridge circuit. Ther- 
mistors have a large negative temperature coefficient 
of resistance, thus providing high sensitivity to small 
temperature changes in the cylinders. Absolute 


' Mev and kev are energy units referring to ‘‘million electron volts” and “‘thou- 
sand electron volts’’, respectively. 








calibrations were made by determining the tempera- 
ture rise when a measured quantity of electric energy 
was dissipated in a resistor embedded in the evlinder. 


3. Construction of the Calorimeter 


Figure 1 shows the vacuum chamber containing 
the two lead cylinders. The cylinders appearing in 
the background are reflections in the chrome lining 
of the chamber. The vacuum chamber was sup- 
ported in a temperature-regulated water bath. 
Figure 2 is a schematic front view showing the posi- 
tion of the vacuum chamber m the water bath. <A 
detailed description of the various components 
follows. 


3.1. Cylinders 


The absorbing cylinders were made of lead because 
the absorption per unit heat capacity for high-energy 


X-rays is greater for a material of high atomic num- | 


ber than for one of low atomic number. The eyl- 
inders were 7.5 em long and 4 em in diameter. Fig- 
ures 1 and 2 show their location in the evacuated 
chamber. The surface of each cylinder was gold 
plated to limit radiation losses, and the cylinders 
were suspended by nylon threads to limit conduction 
losses. The cylinders were cut in planes through 
their axes to simplify the insertion of the thermistors 
and calibrating resistors. Figure 3 shows one of the 
evlinders partly opened. The end of the eylinder 
that received the irradiation was called the front 
surface. When X-rays are absorbed by lead the 
maximum density of energy dissipation is at a posi- 
tion a few millimeters behind the front surface. 
This position depends upon the energy of the incident 
X-rays. The resistors for calibrating were, therefore, 
placed near the front surface. The position of these 
resistors was not critical because the heat from them 
was transmitted throughout the lead eylinder in a 
time that was short compared with the irradiation 
time. The resistors were located near the end of the 
evlinder in order to be out of the X-ray beam where 
they would not affect the absorption of the X-rays 
by the lead. One resistor was located in each half 
of the cylinder. Commercial 1,000-ohm '.-w carbon 
resistors were used. The resistors were coated with 
a liquid cement, placed in a groove in the lead, and 
surrounded with Wood’s metal. The cement pro- 
vided electrical insulation, and the Wood's metal 
filled the air gaps, thus providing better thermal! 
contact. 

The thermistors were located near the back sur- 
face and near the edge of the evlinder, so that they 
also were not in the primary X-ray beam. One 
thermistor was located in each half of the evlinder. 
They were embedded in the lead in the same manner 
as the resistors. 

Number 43 enameled copper wire was used to 
connect the thermistors and the calibrating resistors 
to the outside. The wires were about 10 em long, 
and small enough to limit conduction losses to a low 
value, yet large enough that their resistances were 
small compared with those of the thermistors and 
calibrating resistors. 
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Top view of the twin lead-cylinder calorimeter fo 


Figur: 1. 
measil ing \ -ray-be am power 7 
The lead cylinders in the background are reflections in the chromium-plate 
surface 


3.2. Vacuum Chamber 


The lead evlinders were suspended inside a brass 
vacuum chamber by means of 16 nylon threads 
The chamber had ‘.-in.-thick walls and = '4-in.-thiek 
bottom and top plates. Its inside dimensions were 
6 by 11 in. and 2's in. high. The inside surface was 
lined with 'ye-in. plates, the surface of each of which 
had been chromium plated to limit thermal radiation 
effects on the evlinders. A chromium-plated metal 
partition separated the two evlinders. The chamber 
was evacuated to a pressure of about 107° mm Hg 
by means of an oil-diffusion pump and a mechanical 
forepump. A 1'-in.-diam hole was located just in 
front of each eylinder to permit the X-ray beam to 
enter without penetrating the '-in. wall of brass. 
These holes were sealed with 0.003-in. brass foil as 
its absorption for high-energy X-rays is negligible 
The electric connections were made through a glass- 
to-metal seal containing eight small metal tubes. 
The wires were threaded through but insulated from 
the tubes and then sealed with a liquid cement. 


3.3. Water Bath 


<s shown in figure 2, the vacuum chamber was 
immersed in a temperature-controlled water bath. 
The water container, 15 by 21 by 12 in. high, was 
made of Bakelite reinforced with angle iron to pre- 
vent warping due to water absorption. There was 
at least 3 in. of water on all sides of the vacuum 
chamber. Two stirrers, placed diagonally opposite 
one another, circulated the water. The stirrers 
were geared to s\ nchronous motors in order to help 
maintain 2 constant temperature — distribution 
throughout the water bath. Air-filled Bakelite tubes 
extended through the water bath in front of the 
entrance ports of the vacuum chamber in order to 
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FIGURE 2. Schematic diaeram of a vertical cross section of the calorimeter with the associated water bath 


and equipment, 


keep the water from absorbing energy from the 
X-ray beam. A 0.001-in. aluminum foil was placed 
at the external end of the tubes to limit exchange of 
heat either by convection of air or by radiation to 
the inner foil. 

The vacuum chamber was supported inside the 
water bath on the points of three Lucite cones, each 
of which rested on three marbles. The vacuum con- 
nection to the top of the chamber consisted of a 2-in.- 
diam glass tube on top of a metal sylphon soldered 
to the vacuum chamber. The purpose of this con- 
struction was to provide high insulation between the 
chamber and the outside, and to provide good 
thermal contact between the chamber and the water 
bath, so that outside temperature fluctuations would 
have as small an effect as possible on the temperature 
of the vacuum chamber. 

The temperature of the water bath was regulated 
by means of an immersed helical bimetallic strip with 
an electric contact at one end. The electric current 
through the contact actuated an electric relay that 
controlled the ‘fon’ and ‘off’ time of the electric 
heater at the bottom of the water bath. Thermal 
regulation under optimum conditions was estimated 
to be 0.005 deg C. The heater was operated at 
24 w and maintained the temperature at 33.4° C, 

It was found necessary to surround the Bakelite 
box with a larger plywood box to prevent temper- 
ature fluctuations due to convection currents in the 
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Figure 3. Schematic diagram of one of the lead cylinders split 
in the plane along its axis to show the location of the thermis- 
tors and the heating resistors used for calibration, 


The X-ray beam was incident axially on the front surface. 


room. The temperature of the air in the plywood 
box was controlled by a thermostat to within 40.3 
deg C of 32.4° C. 

The entire calorimeter was mounted on an adjust- 
able hydraulic lift table for convenience in lining up 
the calorimeter with the beam from the X-ray 


souree. 
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Ficure 4. 
a function of temperature. 
rhe temperature coefficient of resistance was —0.0363 for all four thermistors. 


4. Characteristics of the Thermistors 


The temperature-sensing elements in the lead 
cylinders were Veco 3351 thermistors, about \ in. 
in diam by '%. in. long. They were made of a 
mixture of the oxides of cobalt, manganese, and 
nickel. Figure 4 shows a plot of their resistance 
versus temperature over the range of about 30° to 
45° C. Over this range the logarithm of the resist- 
ance is approximately a linear function of the 
temperature. Over a large range the resistance, R, 
at Kelvin temperature 7’, can be expressed as 


R R ea T—1/T,)) 
0 ; 
where Ap is the resistance at temperature 7) (304° K), 


and B is a constant of the material. Table 1 shows 
the characteristics of the four thermistors used. 


TABLE 1. Thermistor characteristics 
rhermistor Po B 
Ohms kK 
l 2,479 3, 355 
: 2,719 3, 359 
3 2, 776 4, 362 
1. 2, 420 3, 353 


5. Wheatstone-Bridge Circuit 


Two thermistors were embedded in each lead 
cylinder and connected in a Wheatstone bridge, as 
shown in figure 5. The 24-ohm copper resistor 
compensated for the adjustable resistance of the 
resistance box and provided an adequate range. 
A mercury cell supplying 1.34 v was used because 
of its stable output voltage with time. A variation 
in the cell voltage would affect the currents in the 
thermistors. These currents affect the tempera- 
tures of the thermistors, which in turn affect their 
resistances, resulting in a possible unbalance of the 
Wheatstone bridge. This effect was small because 
of the similarity of the thermistors, the symmetry 
of the bridge, and the magnitude of the bridge 
current used. The compensating copper resistor 
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Resistance of thermistors used in the calorimeter as 
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Figure 5. Wheatstone-bridge-circuit diagram showing location 


of electric components. 
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and the mercury cell were located in a glass tube 
immersed in the water bath. The tube was stuffed 


with cotton to inhibit temperature variations due 


to convection from above. A 0.01 deg C change i 
temperature of the 24-ohm copper resistor could 
create the same bridge unbalance as a 5X 107° deg C€ 
change in the thermistors in one of the lead cylinders 

The resistance box was a six-decade box with steps 
of 0.01 ohm. Because of unpredictable variations 
in contact resistances and emf’s originating in_ the 
switch, the decade was modified as shown in figur 
6, involving a variation of several methods described 
by Mueller and Wenner [4]. As the contact was 
moved from one end of the decade to the other, the 
total resistance changed from 1.00 to 1.10 ohms 
Because the variable contact resistance was always 
in series with a relatively large resistor, it contributed 
less variation to the over-all resistance. 

The four thermistors, selected from a group of 
12, had the most similar characteristics. This was 
desirable in order to reduce drifts and fluctuations 
at the output of the Wheatstone bridge produced by 
ambient temperature fluctuations and changes in 
bridge working currents that affected the tempera- 
ture and hence the resistance of the thermistors. 
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For : ; 
these expermments, the first two terms of an expan- 
« exponential formula were sufficient, so 
that R could be expressed as 


sion of 


R=R[1— (Bf T— 7To})/(To?)], 

resistance could be assumed to be a linear 
function of the temperature. The temperature 
coefficient of resistance was —0.0363/° C for the 
four thermistors used. As the temperature rise 
during any set of runs did not exceed 0.01 deg C, 
the error introduced by assuming linearity and 
neglecting the next term in the expansion was less 
than 1 part in 10 million. 

Because of the small voltages detected in the 
Wheatstone bridge (order of microvolts), it is neces- 
sary to take every possible precaution to avoid 
variable voltages due to thermoelectric effects when 
conductors of different materials are used. Copper 
connecting wire was used throughout, and a special 
solder with small thermal electromotive force with 
respect copper was used to minimize thermo- 
electric effects. Difficulty was originally observed 
with the use of solid metal leads through the glass 
seal to the vacuum. Although copper wires were 
attached to both ends, the variable thermal electro- 
motive forces due to temperature differences between 
the inside and outside of the vacuum chamber were 
too great to work with. For this reason metal tubes 
were substituted, and the copper connecting wires 
were run through the axes and insulated from the 
metal to form a circuit of copper alone. 

About 0.26 ma of current flowing through each 
arm of the Wheatstone bridge supplied about 0.35 
mw of power continuously to each lead eylinder. 
This caused the equilibrium temperature of the 
cylinders to be of the order of 0.03 deg C above the 
surrounding temperature. As long as this tempera- 
ture difference was constant and all heat losses were 
approximately linear functions of the temperature 
difference, this effect could be disregarded. 

The detector employed in the bridge circuit was a 
d-c breaker amplifier. Its input impedance was 
1,000 ohms. In the breaker amplifier the d-c input 
voltage Was converted to an 8-cps signal by a mechan- 
ical breaker contact. The amplifier, tuned to 8 eps, 
amplified the signal, which was then converted back 
to direct current by a second breaker synchronous 
with the first. In this way the background noise 
was kept low. The observed noise level (maximum 
fluctuations peak to peak) corresponded to about 
0.l-uv input, and the voltage gain was of the order 
of 500,000. 

The temperature data were taken by either of the 
following methods: (1) The output of the amplifier 
was connected to a recording milliammeter. The 
amplifier and recorder were used as a null detector, 
so that variations in gain in the amplifier were not 
of primary importance. As the thermistors changed 
resistance due to a temperature change of the lead 
evlinder, the resistance box was adjusted, as required, 
to bring the bridge back into balance. The required 
balancing resistance was then plotted against time 


and thi 


to 


to provide the data. (2) The second method 
consisted in connecting the output of the amplifies 
to a recording d-c potentiometer (0 to 30 my) wit! 
an 11-in.-wide chart and recording the unbalance ot 
the Wheatstone bridge without changing the balane- 
ing resistor. Although this system eliminated the 
transient voltages caused by switching the balancing 
resistor, it had the disadvantage of depending upon 
the constancy of the amplifier gain. Periodic checks 
were made of the amplifier gain by adding a l-zv 
signal to the input and observing the deflection; vet 
no evidence of a change in gain was observed. 


6. Operating Characteristics 


The smallest short-time signal detectable was 
limited by the noise level of the output signal. This 
noise had an over-all amplitude of about 0.0008-ohm 
equivalent, and corresponded to about 5 107° °C, 

The smallest sustained input power detectable 
was limited by the long-time drifts in the balance 
point of the bridge. These could have been due to 
drifts in the thermistor resistance, drifts in the 
battery voltage, or changes in the water-bath 
temperature that affected the bridge output due to 
asymmetry of the bridge. The observed drift during 
good operation was less than 0.005 ohm/hr, and 
corresponded to approximately 3 107° °C/hr. The 
temperature rise due to the 100-4w input power 
in one cylinder was approximately 2 107° °C /hr and 
could readily be observed and measured. 


7. Heat Losses 


It was desirable to keep the heat losses small for 
several reasons: (1) it put less stringent requirements 
on the constancy of the water-bath temperature and 
led to a greater stability of the bridge balance, (2) a 
larger over-all temperature rise could be accumu- 
lated if the time constant were longer, and (3) linear 
approximations for heat-loss calculations could be 
used. The associated disadvantage was the long 
time required for the calorimeter to attain thermal 
equilibrium. The cooling time constant was de- 
termined from data taken of the required balancing 
resistance versus time before and after the cylinder 
was heated with 14.03 mw for 10 min. The data are 
shown in figure 7. Assuming Newton’s law of cooling 
and linearity between FR and 7, the cooling constant 
K was determined in the following manner: 


dR/idt=—K(R—R,) 


K (dR, /dt—dR,,/dt)/(R; R,), 


where dR,/dt and dR,,/dt are the slopes at R, and R,, 
respectively, and FR, is the equilibrium resistance that 
the exponential curve was approaching. A was 
determined at five points on the decay curve and 
found to average about 0.00548 +0.00002 min“, 
when the bath temperature was 33.4° C. This cor- 
responded to about a 3-hr time constant. As is 
shown below, most of the heat losses were due to 
radiation. Therefore, this decay constant should 
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FIGURE 7 Data from which the cooling constant was calculated. 


vary directly as the cube of the absolute temperature 
of the water bath. Four measurements of the decay 
constant over a temperature range of 5 deg agreed 
with the cube-law dependence on the absolute 
temperature, 

The heat capacity of each evlinder was calculated 
from its dimensions and material to be about 136 
j/°C. The heat loss per unit time per degree tem- 
perature difference between the evlinder and its 
surroundings is equal to the product of the cooling 
constant, A, and the heat capacity. Using these 
values, one obtains 12.8 mw/°C as the rate of heat 
loss by the evlinder to its surroundings. This heat 
loss is believed to be due, primarily, to radiation. 
Because of the roughness of the gold surface on the 
evlinders, the emissivity was not known, but must 
be considerably larger than 0.02, which is the emis- 
sivity for a polished gold surface. 

The rate of heat loss per degree Celsius due to 
conduction was calculated and found to be small 
compared with the observed rate of heat loss. The 
loss due to conduction along the wires was 0.065 
mw/°C, which is only about 0.5 percent of the 
observed value. 

The rate of heat loss per degree Celsius by con- 
duction through the air was 0.017 mw/°C, which is 
only about 0.1 percent of the observed value. The 


pressure used in this calculation was 10™mm He, as 


roe 
ruce ° 


measured with an ion 





Assuming the entire observed rate of heat ‘oss isp 
due to radiation, the effective emissivity of the leag 
eylinder in the vacuum chamber can be calculategy 
from the formula . 

F 

dQ/dt cAe,,;(T'— T} 

where ec is. the  Stefan-Boltzmann — constan 
[5.673 X10-" w/em? per (deg K)*], A is the surfag 
area of the evlinder (120 em®*), ey, is the effectiyy 
emissivity of the evlinder in the vacuum chamber} 
and 7’ and 7) are the absolute temperatures of thet 
evlinder and surroundings, respectively. Assuming 
(T— T)<T) and solving for e,,,, one obtains 


; I dQ 
“ff 4¢eATXT—T)’ * dt 


Substituting the measured value of 12.8 mw/°C for! 
(1/i T— Ty) XdQ/dt), 303° K for To. and the above 
values for ¢ and AA, one obtains a value of 0.17 for thy 
effective emissivity. This is a reasonable value 
considering the matted appearance of the gold sur.J 
face of the lead cylinder. It also indicates that af 
considerable improvement is possible if the gold] 
surface were highly polished. } 


8. Calibration 


The method used for calibrating the calorimeter 
Was to supply a measured quantity of electric energy 
to the resistors embedded in the lead evlinder and 
to determine the resulting rise in temperature of the 
evlinder. The rate of heating and total temperature! 
rise were chosen so as to duplicate as nearly as 
possible the conditions during the measurement 0! 
the N-ray-beam energy. This procedure partial 
eliminates some of the systematic errors in the calor: 
metric measurements. A potentiometer with ant 
accuracy of about 0.1 percent was used to measur 
the voltage across the heating resistor, and a mode 
622 Weston d-c milliammeter with an accuracy o) 
0.5 percent of full-scale deflection was used to meas. 
ure the current. The measurements involved about 
half-seale deflections, so the accuracy of the current 
reading was about 1 percent. 

Figure 8 is a freehand sketch showing the output 
voltage of the Wheatstone bridge when method 1 
described in section 5, was used. The vertical breaks 
represent manual changes of the resistance box of 
0.0L ohm. The values between breaks indicate the 
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Figure 8 Sketch showing a way of recording data by using method 1, described in section 5 A 
The vertical breaks were caused by manual changes in the resistance box of 0.01 ohm. The resistance values are NY 


shown between breaks. 
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resistance: of the box at that time. A change in 
slope of the output versus time in figure 8 follows 
any chance in the power being applied to the cylinder. 

The :sistance required for exact balance could be 


| at any instant by interpolation. This 

was calculated at the instant before each 
break in order to minimize the effect of 
transients that sometimes appeared after a change 
of resistance. Figure 9 shows a plot of the calcu- 
lated balancing resistance as a function of time for a 
1.400-key X-ray run. Straight lines were drawn to 
the average data and to simplify the 
Linearity between the balancing resist- 
ance and temperature was assumed, and the observed 
temperature rise was calc ulated to be 4.0 107* °C 
10 is a graph of a calibration run, using 
method 2 described in see sien 5; 122.6 ww of power 
was applied for 12 min. Linearity between the 
bridge unbalance and the temperature was assumed, 
and the observed temperature rise during the heating 
period in this example was alculated to be 
5.47 X107* °C. 
In figures 9 and 10 there is apparent an approxi- 
mately 45-sec time lag between the application of 
heat to the evlinder and the detected temperature 
rise. This time lag was caused by a combination of 
two phenomena: (1) the time required for apices 
of heat through the lead (most of the heat energy 
was supplied near the front surface, be see the 
thermistors were located near the back surface), and 
(2) the thermal time constant of the thermistor 
determined by its heat capacity and its insulation 
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from the surrounding lead. The correction of the 
observed temperature rise for heat loss was made 


on the basis of Newton’s law of cooling, using the 
observed curve of temperature as a function of time. 
This procedure makes the assumption that the 
average temperature throughout the evlinder is 
equal to the surface temperature, and that the time 
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TABLE 2. Electrical calibration of calorimenter 


Resistance Resistance Resistance 


change cor- change cor- | change cor- 7 ,  Calibra- 
responding responding responding en —- tion 
to tempera- to heat to total pews — factor ! 
ture rise loss energy 
Ohm Ohm Ohm wu min Ohm/pw- 
min 
0. 058 0. 052 0.110 70.6 20 77.6x10°% 
il4 O08 122 69.9 21 83. 2 
. 095 028 . 123 70. 5 22 79.4 
. O80 033 118 70.5 20 80.3 
. O80 . 036 - 116 70.7 20 81.9 
Average =80.5 X10%+1.0X10~ ohm/suw-min. 
lag in the observed temperature is constant. It is 
recognized that these assumptions may not be 


strictly true when comparing the heating portion of 
the curve with the initial- and final-cooling portions. 
However, it is believed that the error arising from 
this assumption is small in comparison with the 
estimated over-all error, and will partially cancel out 
when comparing X-ray energy with electric energy. 

The total energy absorbed by the lead cylinder 
can be expressed as a sum of two terms, one involving 
the observed temperature rise and the other involv- 
ing the energy lost to the surroundings as a result of 
the temperature difference. The heat lost to the 
surroundings was calculated from the average of the 
observed cooling slopes before and after the heating 
period. Table 2 summarizes the results of five cali- 
bration runs over a period of 2 days, using method 1. 
The heat loss was computed for each run, so it was 
not necessary to wait each time for complete tem- 
perature equilibrium to be attained. It is seen that 
at times the heat lost through radiation was just 
about as large as the heat corresponding to the net 
temperature rise. 
this calorimeter 
an rms error of 1.0% 


The average calibration factor of 
with 


was 80.5 107° ohm/uw-min, 
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9. Considerations in Measuring X-ray-Beam 
Power 


The described calorimeter was designed for measur- 
ing X-ray beams in the range 1 to 180 Mev. The 
total beam power will be the sum of the power 
absorbed by the lead cylinder plus that which escapes 
by transmission, scattering, or by other means. As 
the escaping energy must be measured by auxiliary 
experiments, using ionization chambers or scintilla- 


tion detectors (the very methods whose uncertainties | 


make the calorimeter necessary ), it is desirable to 
select X-ray-beam and lead-cylinder dimensions that 


limit the escaping energy to less than 10 percent. If | 


an effort is made to reduce the escaping energy 
further by increasing the dimensions of the lead 
evlinder, the heat capacity would be increased, and 
thus the sensitivity would be decreased. To obtain 
maximum accuracy and sensitivity, a specific size 
and shape of lead cylinder should be chosen for each 
application, depending upon the beam power avail- 
able and the X-ray energy of the source. 


The authors express their gratitude to H. W. Koch) | 
for his encouragement during the development of — 
the calorimeter. The contributions of Samira \M 
Aly * and Donald Bowers during the initial siages 
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Pyrochlore-Type Compounds Containing Double 


Oxides of Trivalent and Tetravalent Ions’ 
Robert S. Roth 


A study 
oxides. 


has been made by 


X-ray diffraction analyses of the A,O;:2BO,-type double 
It was found that many of these mixed oxides, after appropriate heat treatment, 
formed compounds of the formula type A,B,O; 


Most of these compounds crystallized in 


the cubic system with a face-centered cell similar to that found for the mineral pyrochlore, 
although some were distorted from the ideal cubic structure. 


Indexed 
Sim.Oy-2 TiO», 


X-ray 
Cid,O3-2 TiOd, 


diffraction powder 
Dy,Os-2TiO», 


patterns 
Y,05-2Ti0», 


cubic 
La,O,-28nQOnz, 


given for the 
Yb.O3-2 TiO,, 


are compounds 


NdoQOs: 


28nO0,, La,Oy-2ZrO,, Nd.Oy-2ZrO., and for the possible compounds Y ,O 3-2ZrO, and Nd, Os, 


2U0,. 


Unindexed patterns are given for LasO,-2TiOs, Nd,O,-2TiOs and Bi,Os-2SnQb. 


On the basis of the existence of the two compounds La,O,-2ZrO,. and Nd,O;-2ZrOs,, the 
phase diagrams for the systems La,O,-ZrO, and Nd,.O,-ZrO, have been revised. 


1. Introduction 


A partial survey of the reactions occurring in 
binary oxide mixtures of the type A,O;:2BO, has 
heen conducted as part of a program of fundamental 
research on ceramic materials. Combinations of 
simple oxides in this proportion were selected be- 
cause of the current interest in ferroelectric ceramics. 
Cd,Nb,O;, at room temperature, has the cubic pyro- 
chlore structure, although the ferroelectric form at 


low temperatures is distorted from this ideal struc- | 


ture [1, 2, 3]° The fact that the compound 
Cd,Nb.O; is known to be ferroelectric at certain 


| CeO,. Caleined cerium dioxide of nominal 99-per- 
cent purity. 

UO,. Urania, supplied by the United States Atomic 
Energy Commission, of over 99.9-percent 
purity. 

La.O,. Lanthanum oxide of 98-percent purity, re- 

mainder mostly water and CQ,. 

Nd,O;. Neodymium oxide of over 99-pereent pu- 

rity. 

Sm,O;. Samarium oxide of nominal 99-percent 

purity. 

Gd,O;. Gadolinium oxide of nominal 98-percent 

purity. 

Bi,0;. Bismuth trioxide of over 99-percent purity. 


temperatures [1] led to the study of other compounds 


with similar structure. 

At the time this work was initiated, no references 
were known deseribing a pyrochlore structure that 
did not contain a pentavalent ion as an essential 
element of the compound. A ternary compound of 
this type having the formula CaQ-ZrO,-2TiO, or 
(Ca; oZrio) Ti,O; was, however, deseribed by Cough- 
anouret al. [4]. Another reference has recently been 
made by Padrow and Schusterius [5] to pyrochlore 
phases in the system La,Q,-SnO,-TiQ,. 


This report is concerned with some combinations 


of the oxides of trivalent and tetravalent ions which, 
on the basis of radius ratio, might be expected to 
form compounds of the pyrochlore type. 


2. Materials and Methods 
2.1. Materials 


The following substances were used as sources of 
the component oxides in this study: 


MO, A rutile of 99.5-pereent purity. 

SnQ,. Precipitated tin oxide of over 98.5-percent 
purity. 

ZrO... Dense zirconia of over 99-percent purity, 
caleined at 1,450° C. 

' This work has been sponsored as part of a program for Improvement of Piezo- 
electric Ceramics by the Office of Ordnance Research of the Department of the 
Army 

Figures in brackets indicate the literature references at the end of this paper 
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BizOy35n0.-5H,0. Precipitated hydrous bismuth 
stannate of at least 98-percent 


purity. 


Y,O;. Yttrium oxide of 98- to 99-percent purity, 
the remainder being rare earths. 

Dy,O;. Dysprosium oxide of nominal 98-percent 
purity. 

Yb,O,. Ytterbium oxide of nominal 98-percent 
purity. 

InoQO,. Rather impure indium oxide specified “for 
manufacturing use only.” 

Sb.O;. Antimony trioxide of over 99-percent purity, 


with 0.2 pereent of As,Os. 
2.2. Sample Preparation 


The starting materials, in sufficient quantities to 
give either a 10.0-g sample or a 1.0-g sample, depend- 
ing upon the availability of the raw materials, were 
weighed to the nearest milligram. No corrections 
were made for the percentage purity of the raw 
materials, except for loss due to water and carbon 
dioxide. They were then mixed together with a few 
drops of distilled water to assure bonding and formed 
into l-in. or '-in.-diam disks at a pressure of 5,000 
lb/in.*. The pressed disks were fired for 4 hr at 
1,100° © on platinum foil in an air atmosphere, 
using an electrically heated furnace wound with 
80% Pt-20% Rh wire. 















Following the preliminary heat treatment, the 
disks were ground, using a mullite mortar and pestle, 
remixed W ith a few drops of distilled water, and new 
disks, about \ in. high, were formed at 15,000 
lb/in? in either a \-in. or '\-in.-diam mold. These 
specimens were then ready to be used for solid-state 
reaction studies. Some of the specimens used were 
prepared by G. R. Shelton of the Bureau during the 
last 10 vears by methods similar to those described 
here. Most of these specimens were reheated, in 
the present study, to obtain better erystallinity. 


2.3. Study of Solid-State Reactions 


In the study of solid-state reactions, a conventional 
platinum-wound quench furnace was used. The 
temperature in the furnace was measured with a 
Pt versus Pt-10-perecent-Rh thermocouple and con- 
trolled by a modified Roberts-type controller. The 
quenching technique was used because it has been 
observed that sharper X-ray patterns are often 
obtained by very fast cooling of the specimen. 
Phase transitions in the pvrochlore structure are 


probably completely reversible and cannot be 
“frozen in” by quenching. For quench tests, the 


\- or \j-in.-diam disks were placed on a platinum 
platform, which was then suspended in the furnace. 
The test temperature ranged from 1,250° to 1,550° C 
and was maintained constant for a given length of 
time. Equilibrium conditions were usually reached 
in less than 3 hr. Equilibrium was believed to have 
been reached when X-ray patterns of the specimen 
showed only a single phase, or when the pattern of 
a two-phase specimen showed no change with 
successive heat treatment. After this reaction time, 
the sample was quenched in air by lifting the speci- 
men disk on its platinum platform out of the hot 
furnace. The specimen cooled to room temperature 
in less than 2 min. The samples were then examined 
by X-ray diffraction, using a Geiger counter 
diffractometer employing nickel-filtered copper 
radiation. 

In the case of specimens containing UQO.,, an argon 
atmosphere was used and the specimens were not 
quenched. The heating element of the furnace used 
was a molybdenum-wound stabilized ZrO, tube and 
is similar to the one deseribed by Davenport et al. |6}. 


3. The Pyrochlore-Structure Type 


The structure of the mineral pyrochlore and of 
related mineral types has been described by several 
workers [3, 7, 8S]. Pyrochlore has the formula 


NaCaNbTaQ,F, in which Na and Ca can be con- 
sidered the A ions, and Nb and Ta the B ions of 
an A,B.X,; structure. The space group for this 


eubie compound has been found to be Fd3m-—O,,, 
with Z=S8 and a unit cell of 10.35 to 10.40 A. 

In the following description of the pyrochlore 
structure, the cell origin has been taken at a center 
of symmetry and the sets of positions referred to 
may be found on page 341 of reference [9]. In this 
structure type there are 16 A ions in position (¢), 
16 B ions in position (d), and 8 negative tons in 
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position (a); for the mineral pyrochlore thes: are f 
the fluorine ions. 
ions (48 oxygen) are in position (f), which contains f 
only one unknown parameter. From a consideration 
of the spatial requirements in the pyrochlore si rue. 
ture, it is obvious that the value of the unknown. 
parameter must be equal to or very close to *. Tf 
the value of * is used, then all possible ryz combi. 
nations of 's, %, %, and 4, are occupied by oxygen, 
except the eight sites required by the special position 
(b), which are vacant. In such a structure, if these 
vacant sites were also occupied, the formula type 
would be A,B,O, or ABO, instead of A,B.O;. The 
pyrochlore-structure type may therefore be com. 
pared to a fluorite structure with 1 out of every § 
negative ions missing. From data given by Bystrom 
[3], the value of x in the oxygen position for Cd,Nb,0, 
and Cd,Ta,O; can be taken as 0.414 instead of 0.375 
for *%. It will be seen in the following discussions 
that the value of x is probably not the same for all of 
the pyrochlore-type compounds. 

In addition to the extinctions to be expected from 
a face-centered cell, there are two special extinction 
rules due to the particular positions of the ions in the 
pyvrochlore-tyvpe structure: 

(1) The positions (a) and (f), filled by negative 
ions, are such that possible reflections are limited to 
the tvpe hk/ where h+-k+/—2n +1 or 4n. 

(2) The positions (c¢) and (d), filled by positive 
ions, are such that possible reflections are limited to 
planes of the type 


hkl ch 2n- l tn —--2 a7 








y, J 
The result is that, in the pyrochlore-tvpe diffraction 
patterns, certain reflections are always absent. Of 
those that may be observed, some are due only to 
metal atom seattering, some only to oxygen atom 
scattering, and the remainder to both metal and 
oxygen seattering. Therefore, the intensities of the 
peaks in the X-ray powder diffraction patterns are 
characteristic of the pyrochlore structure, and are 
evidence for the existence of such a structure. — In- 
tensity calculations have been made for some of the 
cubic compounds studied and will be discussed in 
later sections. 


4. Results and Discussion 


4.1. General 


The X-ray diffraction powder patterns prepared in 
this study revealed that many of the mixtures 
formed binary compounds of the tvpe A.B.O;. The 
great majority of the compounds were cubic, essen- 
tially isostructural with the mineral pyrochlore. 
Table 1 lists the indexed X-ray diffraction powder 
patterns for these cubic compounds, All of the hél 
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Heat treatment and structure type of compositions studied 


Heat treatment 


Temper- 


ature 
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5a 
5M) 


$25 
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5A) 


WO) 
20 
5M) 
m4) 
nna) 


= a4) 
30 


600 
rn 


Lime 


Structure type 


Distorted pyrochlore 
ao 


Unknown 


ao 


Symmetry 


Pyrochlore Cubie 10, 22s 
do do 10. Is] 
do do 10. 1 
do do 10. 093 
do do 10. 030 
No compound 
do 
Pyrochlore Cubic 10. 702 
do do 10. 568 
Distorted pyrochlore Unknown 
Pyrochlore Cubie 10. 793 
do do 10. 648 
do do 10. 402(? 
No Compound 
do , 
Fluorite solid solution (?) Cubic 449 


No compound 


* Radius of the ions taken from Green [10] 
All specimens, unless otherwise stated, were quenched from the des 
hese specimens were heated in an argon atmosphere to pre 


values allowed by a face-centered lattice have been | diagram) 


listed, up to and including A?+-4?+-/?=192. It may be 
seen that none of the peaks forbidden by the struc- 
tural arrangement of the pyrochlore type are ob- 
served in the X-ray patterns of any of the compounds 
studied. 

Table 2 lists the lattice parameters of these com- 
pounds, together with the composition, heat treat- 
ment, and structure identification for all the materials 
examined, All compounds and mixtures studied 
have been plotted in figure 1 on the basis of con- 
stituent ionic radii. The values for the radii of the 
ions have been taken from Green’s Geochemical 
Table of the Elements [10]. The radius of the A* 
ion is plotted as the ordinate and that of the B** ion 
as the abscissa. It can be seen that only the larger 
A ions with the smaller B ions (the upper left of the 
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to the constituent ionic radit. 


Pyrochlore (unknown symmetry), 
Pyrochlore (cubic). 


Py Single-phase face-centered cubic. 
No compound, 


vent the oxidation o 
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20 


gnated temperature 


f UOo and were t quenched 


form compounds, In the lower-right 
portion of the diagram no compounds are indicated 
Thus, areas of compound formation and also of 
symmetry types can be located on this diagram 
This is similar to the types of diagrams given by 
Wood [11] and by Keith and Roy [12] for perovskit 
structures. Detailed discussion of each of the com. 
pounds studied is given in the following sections. 


4.2. Titanates 


Most of the rare-earth titanates form compounds 
having the cubic pyrochlore structure; however 
La.Ti.O; and Nd.TiO; are not cubic. The X-ray 
diffraction powder data for these two compounds 
are listed in table 3. The X-ray pattern for 
La, Ti,O; is very sharp and clear; however, that of 
Nd, Tif ), is rather diffuse. These compounds seem 
to be approximately isostructural with CaZrTi,0. 
4]. The X-ray patterns of these compounds al 
have diffraction peaks more or less grouped around 
the positions that would be expected for a cubi 
pyrochlore structure. An effort has been made t 
correlate the structure of these titanates with the 
orthorhombic form of the A,B,.X; compounds exem- 
plified by the mineral weberite, as described by 
Bystrom [7] for Ca,Sb,O; and other antimonates 
The titanates reported in the present study do not 
seem to have the same structure as that reported 
for the compounds studied by Bystrom, and the X-ray 
patterns strongly suggest distortions of the pyro 
chlore type. 

The other four rare-earth titanates 
have the cubie p\ rochlore structure, as does vttriun 
titanate (table 1). These compounds all give very 
sharp X-ray patterns, with relatively strong addi 
tional lines for the Akl peaks required by the pyro 
chlore structure as compared to the fluorite structure 
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TABLE 3. X-ray diffraction powder duta for the comp 
La,O;-2TiO, and Nd,O;-2TiO.— Continued 


LayOs-2TiO; Nd)Oy 2Ti0O2 
29 d ]* 20 d I 
1.2 1. O7S 
) 1. 07 7 
uw? Mw) 1. 066 Q 04.06 1.053 ( 
U3. OY 1. O56 13 95. AT 1.040 4 
Ww. 72 1.031 ' 
05.13 1.044 5 47. 20 1. 027 
W686. 15 1. 035 7 97. 82 1. 022 t 
US. O04 1. 020 4 
YS. US 1.013 3 100. 70 1. 000 
100. 45 1. OO2 6 i] > 2 , 
100. 83 0. 9905 e { 103. 12 0. OS34 
P eyapen f 106.02 9643 i 
103. 96 9 1 106. 60 O607 7 
104. S88 V717 3 | - = as 
105. 22 O64 r { 107. 70 9539 
108, 40 W407 4 110.79 US45S 4 
110. 25 U3SS8 ot) 112.13 W224 
111.93 W204 6 
113.31 9221 7 116. 45 9061 7 
113.84 v192 s 116, 87 O40 2 
119. 03 . SSS 5 
120. 40 . S876 6 
120. 70 SSH 7 124. 02 8725 6 
123.15 S758 i 
124. 0S S720 1 
124. 68 SHO { 
125. 32 S671 1 
128. 35 S557 4 
128. 81 S540 3 
130.14 S404 5 
130, 82 S471 5 
133. 35 SASS S 139. 41 S212 0 
135. 10 8334 5 
135. 80 S313 Ai) 
137. 80 S25) ‘ 149. 87 . 7976 
147.10 SOSI 6 
148. 02 sol2 
149. 92 7976 5 ° 
150. 45 TOH6 s 158. 45 7M t 
151.65 744 7 
152. 00 7923 


» J is the observed height of the diffraction peaks 
hese peaks were too strong to be completely recorded on the X-ray pattern. 
Table 4 shows the intensity values calculated for a 
few representative planes of Sm,7Ti,O; and other 
compounds. The intensity valves were calculated, 
using the formula 
y ] T cos” 26 
Toate. ~# ? sin?6-cosé’ 
where 
F'=the structure amplitude as given in reference 
19, p. 519]. 
p=multiplicity factor, and the balance of the 
equation refers to the polarization and 
Lorenz factors for powder patterns. 
These intensities were then recalculated on the basis 
of 100 for the strongest peak. ‘Table 4 shows the 
individual contribution of the ions in the structure 
and the final intensity valves compared with the 
measured peak heights. It can be seen that both the 
0.375 value and the 0.414 value of the z parameter 
for 48 oxygen give poor agreement between observed 
peak heights and calculated intensity for Sm,Ti,O, 
although the agreement of the zirconates and stan- 
nates is fair for 0.414. It must be concluded that 
the oxygen parameter for the titanate compounds is 
slightly different than that reported by Bystrom [3] 
for the niobates and tantalates. 








FABLE 4 Calculated versus observed intensities for some of the 
i: ffraction peaks of typical A,** B,**O; compounds 
‘x4 O Calculated J * 
16 Observed 
" R iws) I . 
0. 375 0.414 r=0. 375\r=0. 414 
Sm. ThO- 

1! Hie 126 4s 0 70 5.6 33.1 22 
>») R20 242 0 0 0 100.0 100.0 b> 1004 
wn) “ 73. 22 iY nO 2H. 1 35. 1 “wo 
$31 Se 112 30 0 6 0.6 21.7 3b 
11/333 rp) au 52 0 3l 0.9 0.3 14 

Nd)ZrmO 
111 2 73 is 0 7u Ww. 5 46 5 
222 sid 514 0 0 0 100.0 100.0 » 100+ 
TED TSH ri 23 274 153 26. 7 33. 6 is 
$41 isd 245 sl 0 Os :. 2 58 s 
11,333 25 15s 53 0 3l 4.1 3.0 
YoZr0 
11 271 24 is 0 70 0.9 “03 0 
22 yuu 512 a “ 0 moo to 1+ 
Hi) is2 i) 23 274 153 24.6 34.1 i3 
$31 238 246 3l 0 is 1 15 0 
511.333 Hiv i 53 0 31 2 Os “ 
NabsnmO 
lll ie 300 is 0 70 3.3 og { 
222 sil tol 0 0 0 10.0 100.0 100+ 
win) 782 (42 23 J74 153 31.5 ss. 77 
$31 sav 313 31 0 Hs Os 22 2 
511 333 TA SSu 53 0 1 { “02 4 


* J is the height of the diffraction peaks 


ihese peaks were too strong to be completely recorded on the X-ray puttern 


Neither indium nor antimony oxides form com- 
pounds with titania. The X-ray patterns for the 
mixtures In,O,2TiO, and Sb,O,:2TiO, show only the 
presence of TiQ,. The specimens had contracted 
and lost weight, indicating that the other oxide had 
volatilized and had not combined with titania. It 
was found that, if In,O, or Sb.O, reacts with a second 
oxide to vield a true compound, little or no material 
is lost in volatilization under the heating conditions 
used in these experiments. 


4.3. Stannates 


Both of the rare-earth oxides, lanthana and neo- 
dymia, form cubic compounds with stannic oxide. 
There are only a few small peaks in the X-ray dif- 
fraction powder patterns (table 1) that cannot be 
indexed on the basis of a fluorite-type structure with 
a unit-cell size of approximately 5 A. However, 
these additional lines indieate that the true structure 
is that of the pyrochlore type with a cell size of ap- 
proximately 10 A. The diminished intensity of 
these excess lines, as compared to those of the tita- 
nates, is due to the small difference in the seattering 
power of the two metal ions, as shown in table 4 
for Nd.Sn.O;. No other rare-earth stannates have 
been studied, but it is probable that many of them 
form cubie pyrochlore compounds. Also, it is evi- 
dent, from the work of Padrow and Schusterius [5] 
that solid solutions of the rare-earth titanates and 
Stanhates can oecur. 
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TABLE 5. X-ray powder diffraction data for the con ound ; 


Bi,O,.28n0, 





d Ta d Ja d I 
ad 2 

‘ 
6.2 2 1. 496 7 1. 028 ] i 
3.77 14 1. 428 5 1. 006 i 
3. 61 1) 1. 425 in 0. 9964 { ' 
3. 22 27 1.419 10 . 9924 
3. OS » 100+++44 1.304 7 ORl4 
2. 78 2! 1.391 4 O76 t 
2. 67 b 10044 1. 372 ' W641 { : 
2.47 a 1.349 s O446 f 
2. 45 22 1. 336 17 338 ‘ 
2. 39 s 1. 328 t W315 

t 
2. 35 13 1. 206 ‘ W277 
2. 240 10 1. 259 1 9200 { 
2. ISI 5 1. 2453 s G168 } 
2. 090) 5 1. 235 7 WOOT t 
2 O65 7 1. 226 30 W033 M4 
2. 057 11 1. 205 ’ SOS ’ 
1. 448 nH 1.105 2 SAIT i 
1. SSS 100+ 1.175 7 S7H5 H 
1. S68 s 1.173 7 870, ‘ 
1.812 21 1. 166 ‘ S45 

' 
1. 807 1 1. 161 S5N4 
1. 780 1! 1. 148 t S4N0) ‘ 
1. 7#3 14 1. 139 ' S444 s 
1.717 rt 1. 122 8374 * 
1. 690 s 1.110 ‘ S346 : 
1.629 ” 1.091 1 8242 
1.610 row 1.075 S171 s F 
1. 597 12 1.04 i S144 ? 
1. 42 32 1.045 j SO4 rn 
1. 531 13 1. 033 ’ JUST s 


* Jis the observed height of the diffraction peaks 
» These pe iks were too strong to be complete ly recorded on the X-ray pratterr 


The reaction between Bi,Q, and SnO, has been} 


studied extensively by Coffeen [13, 14], who indi 
eated that the hydrous bismuth stannate, Bi, 
(SnQ,),-5H.O, fired at 1,149° C, formed the com-% 


pound Bi,(SnOs)s. The present study indicates that} 
the Bi,O,:3Sn0, mixture contains a new compound 
plus free SnO,, and that the composition Bi,O,- 2Sn0,) 
contains no free SnOQ,. It had been reported by 
Aurivillius [15] that the system Bi,O,-TiO, contained 
a compound, Bi,Ti,0,., that was pseudotetragonal 
actually orthorhombic; therefore, the analogous} 
compound 2Bi,0,:38n0, was looked for in the! 
Bi,O.-SnQO, svstem. The reported structure could 
not be correlated with the compound found in the§ 
Bi,Q,-SnO, system. It must therefore be concluded 
that the Bi,O.-SnO, compound is not isostructural 
with Bi,Ti,0,., but has a distorted pyrochlore-type 
structure with a ratio of Bi,O,:25nQk,. p 

The X-ray pattern of the Bi,O,-28nQO, compound) 
given in table 5 resembles very strongly those of the} 
cubic pyrochlores. All of the cubic peaks for a 
10.68 A cubie pyrochlore structure can be found i 
this pattern, although there are other peaks present! 
A cubie unit cell of 21.37 A could account for mosth 
of these excess peaks, but some would be still unex-/ 
plained. Although the pattern is) very” strongly) 
pseudo-cubic, it must be concluded that the true) 
symmetry is other than cubic. 

It may be seen in figure 1 that the radius of the 
Bi*® ion is close to that of Gd** and, on the basis off 
radius ratio only, could be expected to form a cubicg 
pyrochlore structure with SnO,. However,  thel 
electronic polarization of the Bi® ton is probably 
different from that of the rare earths and may bef 
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This notion of electronic polarization has been very 


related to symmetry. 
valualic in explaining the distortion from cubic 


symmetry in the perovskite structures [16]. 
4.4. Zirconates 


Specimens of La,O;:2ZrO, and Nd,O;:2ZrO, 
were carefully examined by X-ray diffraction powder 
analysis because previous work by Brown and 
Duwez [17] indicated that no compounds existed 
in these two systems. ‘Table 1 shows that several 
extra lines of relatively low intensity were found in 
both patterns, which could be indexed only on the 
basis of a cubic pyrochlore-type compound. In 
table 4 the intensities calculated for Nd.Zr,O,;, using 
an oxygen parameter of 0.414, agree very well with 


the observed peak heights. A compound at a 
composition of La,Q,:2ZrO, was suggested by 


Trombe and Foex [18]. From the data presented 
it must be concluded that the compounds La,Zr,O,; 
and Nd.Zr,O; do indeed exist and probably have a 
wide solid-solution range. They would fall in the 
middle of the stabilized cubic zirconia solid-solution 
regions, shown in the phase diagrams of the systems 
La,O,-ZrO, and Nd.,O,-ZrO., proposed by Brown 
and Duwez [17] and reproduced here as figures 2 (a) 
and 3 (a). These compounds are indicated in the 
proposed revision of the two phase diagrams, 
figures 2 (b) and (ec), and 8 (b). The compounds 
LaOy2ZrO, and Nd,O.-2ZrO, have been in- 
dicated as melting congruently in figures 2 (b) and 
3 (b) for the purpose of simplicity of the phase 
diagram, although no information is available on 


the melting points. If the compound melted in 
congruently, there could be no field of La,Zr,O, (or 
Nd,.Zr,O;) solid solution on the high ZrO, side of the 
compound. This would not agree with the 
served widening of the single-phase cubic solid- 
solution area shown by Brown and Duwez [17]. 
Figure 2 (c) indicates the phase equilibria in the 
La,O;-ZrO, system if the pyrochlore compound 
actually melted incongruently. 

It might be thought possible to have both a com- 
pound and a cubic ZrO, solid solution in the phase 
diagrams. Such a situation would demand that 
the diffraction peaks in the X-ray patterns of the 
fluorite-type solid-solution phases coincided with 
the strong peaks of the pyrochlore compound, so 
as to agree with the observed single-phase solid 
solution. A highly complex diagram of this sort 
could be drawn for either a congruently or an in- 
congruently melting compound. As a complexity 
of such a nature is unlikely, it is thought that the 
phase diagrams, figures 2 (b) and 3 (b), best fit all 
of the observed data. 

It has been claimed by Ruff and Ebert [19], and 
Duwez, Brown, and Odel [20] that Sc,O3, Sm,Os, 
Gd.O;, and Y,O, all form cubie zirconia solid solu- 
tions. These conclusions have been quoted by 
Dietzel and Tober [21]. In the present study, 
only the Y,O,:2ZrO, composition has been exam- 
ined. This specimen is apparently that of a single- 
phase fluorite material with no additional lines 
observable. 

It is possible that the Y,O,:2ZrO, composition 
is a true pyrochlore compound in which the faint 
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FiGurE 2. Possible phase equilibria in lanthana-zirconia system, 


Phase 


showing 


equilibria in the system lanthana-zirconia according to Brown and Duwez [17] 
the compound LajZr2O0; melting congruently with solid solution on both sides of the compound 


(b) Postulated phase equilibria in the system lanthana-zirconia showing 
(c) Possible phase equilibria in the system lanthana-zirconia 


the compound La2Zr2O0; melting incongruently with solid solution only on the high La2Qg side of the compound. 
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Possible phase equilibria in the neodymia-zirconia system. 


FIGURE 3. 


(a) Phase equilibria in the system neodymia-zirconia according to Brown and Duwez [17] 
in the system neodymia-zirconia showing the compound Nd2Zr2O; melting congruently with solid solution on both 


sides of the compound. 


lines characteristic of pyrochlore are too weak to be 
observed. As shown in table 4, the intensity of 
these lines is almost entirely due to the difference in 
the scattering power of the A and B ions. The 
relative intensity of the additional peaks should be 
very small as the scattering factors of Y** and 
Zr** are almost exactly equal and the intensities 
depend almost entirely on the oxygen ions. The 


observed peak heights and the calculated intensities | 


of Y.Zr,O; show good correlation. Although it 
cannot be unequivocally decided whether the 


Y.O,;:2ZrO, composition is a fluorite solid solution 
or a pyrochlore-type compound, the latter is more 
probable, and it is not surprising that the excess 
lines cannot be seen on the X-ray diffraction powder 
patterns. It is, therefore, proposed that no sta- 
bilized cubie ZrO, solid solutions exist in mixtures 
of ZrO, and the larger trivalent ions. All such 
reported occurrences would than be solid solutions 
of A,B.O; pyrochlore compounds, as shown in 
figures 2 (b) and 3 (b). 

The In,Qs:2ZrO, composition also shows only a 
single cubic fluorite-type phase in the X-ray diffrac- 
tion powder pattern. As the In,Q,-ZrO, system is 
known to contain two phases in this area [22], and 
the present specimen contracted considerably and 
lost weight, it is assumed that the excess In.O, 
was volatilized and that a cubic ZrO, solid solution 
was actually found in this case. 


| 
| 
| 
| 


24 





(b) Postulated phase equilibria 





A general rule can now be given for the reaction | 


of ZrO, with various oxides of trivalent metals: 
Solid-state reactions of ZrO, with oxides of the 

smaller trivalent ions, for example In,O,, yield cubic 

ZrO, solid solutions; however, the larger trivalent 





ions probably result in A:B.O; pyrochlore-type 
solid solutions. 


4.5. Cerates and Uranates 
a. Cerates 


The only cerate composition examined in_ this 
study was In,O,:2CeO,. The X-ray pattern showed 
only CeO, with no solid-solution or compound forma: 
tion. The In.O, was believed to have volatilized| 

' 





b. Uranates 


Only two compositions containing UQ,, applicabl 
to the present work, Nd.O,:2UQ, and Y,O,:2U0, 


were examined. These preparations have beer 
described elsewhere [23] and will be reviewe 
briefly here. The Y.O,:2UO0O, composition way 


definitely two phases, showing partial solid solutior 
of both components. The X-ray pattern of th 
Nd,O;:2U0, composition showed a single phase 
apparently a fluorite-type solid solution (table 1 
In this case the characteristic lines of the pyrochlon 
structure should easily have been observed, as Ne 
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has an «tomic number of 60, whereas U is 92. When 
the oxygen positions were shifted considerably from 
those of the other pyrochlores in an attempt to 
account for the zero intensity of the additional 


peaks, improbable uranium to oxygen bond lengths 
resulted. It must be concluded, therefore, that a 


compound of the type Nd,U.,O; probably does not 
exist and that the composition Nd,O,:2UQ, is a 
single-phase fluorite-type solid solution with an 
oxygen deficiency. 


5. Summary 
X-ray diffraction powder data have been presented 


to show that many A,O,;:2BO, compositions yield 
pyrochlore-type compounds when heated to the 


appropriate temperatures to give solid-state reac- | 


tions. Oxides of the largest A*® ions with TiO, 
give distorted pyrochlore structures, whereas most 
of the oxides of the large A** ions form cubic pyro- 
chlore compounds with the tetravalent ions. Oxides 
of the smaller A** ions do not form compounds. 
This information is shown in a plot of the radius 


of the A** ions versus that of the B** ions. Here it 
ean be seen that certain areas of this diagram 


enclose certain structural types of compounds, 
although exceptions may be found. 

Because of the discovery of small peaks character- 
istic of the pyrochlore structure in the X-ray dif- 
fraction powder patterns of La,O,-2ZrO, and 
Nd,O;-2ZrO,, these compositions have been called 
pyrochlore-type compounds. The presence of these 
compounds necessitates a revision of the published 
phase diagrams for these systems, and such revisions 
have been shown. 
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Thermal Decomposition of Polytetrafluoroethylene 


in Various Gaseous Atmospheres ”’ 
L. A. Wall and J. D. Michaelsen 


The effects of various gaseous atmospheres on the thermal decomposition of polytetra- 


fluoroethylene have been investigated in the range 450 
The inhibitory substances, chlorine, hydrogen, carbon tetrachloride, 


tion were observed. 


to 500° C. Both catalysis and inhibi- 


and toluene, all produced marked changes in the mechanism of polytetrafluoroethylene 


decomposition. 


Whereas the normal degradation reaction in a vacuum produces mainly 


monomer, the inhibited reactions give large fragments plus side reactions, such as the de- 


fluorination of the skeletal carbon structure 


The results can be interpreted in terms of 


previous theoretical treatments of chain mechanisms for polymer decomposition and lead 
to the conclusion that, although the inhibited reactions produce little monomer and are 
highly random, the relatively high rates observed indicate induced decomposition by the 


inhibitory substances. 
the postulated mechanisms. 


1. Introduction 


The thermal decomposition of poly tetrafluoro- 
ethylene in vacuum has been studied in some de- 
tail by Madorsky and coworkers [1]°. Although 
Madorsky et al. report that the polymer degrades 
by a first-order rate law, close inspection of their 
data (fig. 1) reveals that this is only strictly true at 
temperatures above about 510° C, whereas below 
about 480° C the polymer degrades by a zero-order 
law. 

The existence of the zero-order law governing the 
kinetics of degradation at temperatures lower than 
about 480° ( has been confirmed by one of the 
present experiments, the data of which are plotted 
in figure 2.) Here, the weight-loss-versus-time curve 
for the degradation performed in a flow system under 
nitrogen is given. The plot shows that after an 
initial slight dip the rate is constant up to about 40 
percent of volatilization at 460° C. In the temper- 
ature range studied, however, the energy of activa- 
tion obtained from extrapolated initial rates was SO 
keal, as reported by Madorsky [1]. 

During the degradation the polymer retains its 
rigid form about 50-percent conversion, after 
which it softens. During this pyrolysis, close to 
100 percent of the volatile matter is monomer. 


to 


Most polymer decompositions can be treated in 
terms of a mechanism analogous to polymerization 
involving four elementary processes: initiation, prop- 
agation, transfer, and termination [2]. In the case 
of this polymer, however, the characteristics of the 
rate versus conversion curves (fig. 1) and the type 
of products obtained from the decomposition lead 
one to the conclusion that a chain mechanism 
operative and that little transfer occurring. 


Is 


Is 
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Energetic considerations are discussed and appear consistent with 


Therefore, only three of these reactions have been 
considered here: 


P—+R—C F.—CF,;-+R,—CF:>- 1) Initiation 
(m 1 or 2) 
R—CF,—CF.:—3R-+ C.F, (2) Propagation 
2R—C F.-— 3P or 2P 
(3) Termination 
H 
or R (] seVvaporation J 


The P represents any size polymer, R is part of a 
polymer chain, and the #’s are the specific rate con- 
stants. 

The mechanism then is one in which a bond breaks, 
probably at the ends, forming radicals that unravel 
into monomer and eventually terminate by one or 
more of the processes listed. 

The amount of unraveling per initiation or ter- 
mination act is the kinetic chain length, Y, and a 
variation in this can influence the order of the rate 
with respect to mass, as shown in figure 3. Here are 
plotted theoretical curves, the upper two showing 
the over-all rate for this type of process, in which a 
polymer had an initial degree of polymerization of 
about 1,000 and a kinetie chain length of 10* and 
10°, respectively, with a given rate of initiation. 
With the same rate of initiation but now restricted 
to end initiation, and with little or no unraveling, 
i. «., Y=1, the dashed curve would apply, which has 
about 0.001 the initial rate of the upper curve. This 
dashed curve gives a zero-order rate with respect to 
mass, which would be the situation if complete 
inhibition could be produced. 

The change in order of rate of degradation of 
polytetrafluoroethvlene from zero to first order with 
increase in temperature may indicate an increase in 
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Figure 1. Rate of thermal degradation of  polytetrafluoro- 


ethylene in a vacuum as a function of percentage of volatiliza- 
tion [1]. 
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TIME. HOURS 


Thermal degradation of polytetrafluoroethylene in 


flow system under nitrogen at 460° C 


Figure 2. 


kinetic chain length from a value much lower than 
the degree of polymerization at temperatures of 
degradation below about 480° C to a value equal to 
or greater than the degree of polymerization at 
temperatures above 510° C. A _ possible alternate 
explanation could be that a surface effect is important 
at the lower temperature range, such that the rate 
is proportional to the surface area of the degrading 
saline. whereas at high temperatures the rate 
eventually becomes proportional to the bulk of the 
material. However, in the experiments under gases 
reported here, surface effect does not appear to exist 
to any great extent because the surface area was 
varied in this study. 

The curve with a maximum in figure 3 represents 
the behavior for a purely random degradation of a 
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Figure 3. Rate of volatilization as a function of conversion, 
and kinetic chain length with random initiation (solid line 
and end initiation (dashed line) {2}. 


polymer. Random degradation may be the result of! 
only random initiation, or of a high number of} 
transfer reactions, if in either case a small number of! 
depropagation or unraveling acts occur. A mudd 
slower over-all rate should result if the random effee 
is due to the initiation step alone. 

Whenever a reaction proceeds by a chain mech-| 
anism it should be possible, in theory, to inhibit th 
over-all reaction by adding a reagent that would inf 
effect make reaction 4 in the above mechanism mor 
probable than is normal for the uninhibited reaction 
Changing the history of the polymer in regard t 
catalyst and structural contaminants has previousl 
been shown to have little effect on the rate of degra. 
dation [3]. However, in the thermal decompositio 
of hydrocarbons, various conditions of inhibitior 
have been studied by Rice and Polly [4] and Hinshel- 
wood and coworkers [5]. For gas-phase reaction 
inhibition has been shown to be produced by som 
e. NO and CH;CHCHs. In _ soli 
polymers, and especially in polytetrafluoroethylene 
one of the difficulties would be poor accessibility o 
the inhibitor to the radical chain carrier. Thi 
difficulty was borne out by the failure to alte 
appreciably the course of its degradation in over 4 
experiments in which solid reagents that were con 
sidered to be potential inhibitors were intimatel 
mixed with the polymer. On the other hand, gase 
can diffuse much more readily throughout a reactio 
medium and hence offered, on a priori considerations 
the best possibility of influencing the course of th 
reaction. 


gases, g., 


2. Experimental Reagents 


The polymer used was in the form of a wate 
suspension obtained from E. I. duPont de Nemours ¢ 
Co., Inc. and designated by them as TD-3. Th 
suspension was frozen and then dried in a vacuum 
and the powder so obtained was employed in th 
majority of the experiments. 

The gases used were the purest obtainable fron 
the Matheson Co., Inc. The toluene was Baker ¢ 
Adamson’s ACS grade. 
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3. Procedure 


The rate of degradation was measured by deter- 
mining the weight loss of the polymer after periodic 
heating under the various flowing gases. The 
sample of polymer was usually about 1g and was 
placed in a porcelain boat for heating. The boat 
was heated in a horizontally placed hard-glass tube 
surrounded by a brass furnace that covered about 
one-third of the tube. The gas and boat were intro- 
duced into one end of the tube, and the gas was 
exhausted at the other end. The boat was inserted 
into the tube via a female 34/45 standard-taper glass 
eround joint, after which the tube was capped with a 
male 34/45 standard-taper ground joint. A heavy 
nichrome Wire passing through a small tube in the 
cap, and sealed with a rubber stopper, enabled the 
operator to move the porcelain boat while in the 
tube in and out of the furnace. 

In a typical measurement the boat containing a 
weighed sample of polymer was placed in the un- 
heated portion of the tube. The system was then 
thoroughly flushed with the gas under study and a 
steady flow maintained. The temperature of the 
furnace Was brought to the desired level. The boat 
was then pushed into the furnace. The manually 
operated variable transformer sometimes had to be 
readjusted in order to maintain the desired temper- 
ature. Temperature equilibrium was usually reached 
about 3 to 4 min after the introduction of the boat 
into the furnace. The time at which the tempera- 
ture reattained the desired value was taken as the 
zero time for that particular run. At the end of a 
given interval of time (usually 1 hr) the boat was 
withdrawn from the furnace into the unheated part 
of the tube. This time was taken as the end of the 
interval of pyrolysis. Experience showed that it 
made no difference whether the boat was then with- 
drawn from the tube and allowed to cool to room 
temperature in the air or left in the tube to cool. 

The temperature of the furnace was measured by 
means of an iron-constantan thermocouple read on 
a potentiometer. The thermocouple was centered 
in the furnace by means of a thermocouple well, 
securely attached to the cap. The temperature was 
determined at a position above the boat and was 
steady to +2 deg C. 

When pyrolyzed in a vacuum or under nitrogen, 
polytetrafluoroethylene retains its original form. In 
the experiments reported here the polymer, with one 
exception, became very fluid at the temperature of 
pyrolysis. Although the molecular weight of the 
polymer cannot be determined quantitatively, this 
change in viscosity from a solid to a very fluid liquid 
was taken as an indication of a very pronounced drop 
in molecular weight during the course of pyrolysis 


under the gases. | 
4. Results 


The studied may be divided into three 
groups: (1) those exhibiting a strong catalytic effect 
on the rate of degradation, i O,, NO, H,S, and 
SO,; (2) those producing an initial inhibitory effect, 
1.e., H,, Cl, CCl, and toluene; and (3) those that 


cr var™ 
rases 
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show neither effect, ie., N, and benzotrifluoride. 
The rates under nitrogen (e.g., 0.017%/min at 457 
C) are quite close to those observed previously in 
vacuum work (0.017%/min at 453° C) [1]. Conse- 
quently, results under nitrogen are used for com- 
parison purposes. 


TABLE 1. Gaseous products from the decomposition of poly- 


tetrafluoroethylene under various gases 


Gaseous atmosphere 


Product 


Or Hy, Toluene Cla 
Mole © Mole © Mole ©; Mole ©; 

Perfluoroethylene 33 45 (4) 
Silicon tetrafluoride 19 34 i) (*) 
Carbon dioxide . 63 33 23 (*) 
Carbon monoxide 18 . (*) 
Trifluoromethane 18 (s 
Perfluorocyclobutane __ _. 5 (®) 


* Not identified. The products, besides an excess of chlorine atoms, showed 
ions in the mass spectrograph containing several chlorine atoms as well as oxygen 
and silicon. For example, the following ions were tentatively identified: CCly, 
SiC], SiICh Fe, OSICCl;, and SiCeClk. The largest of these peaks, mass 161, was 
11.7 percent of the maximum peak (mass 70, that is, Clo"). 


4.1. Catalytic Effect 


The curves representing the degradation of poly- 
tetrafluoroethylene at 457° C under the various 
gases exhibiting 2 catalytic effect are shown in 
figure 4. 

In the degradation under oxygen the polymer re- 
mained white and solid throughout the entire con- 
version. This fact, plus the catalytic effect and the 
absence of monomer in the products lead to the 
belief that the gas attacks the free radicals and in- 
creases, in effect, the rate of propagation. The 
analyses obtained from the mass spectra of the gases 
evolved from this experiment are listed in table 1. 
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Figure 4. Weight loss of polytetrafluoroethylene at 457° C, 


under various gases having a catalytic effect on the thermal 
degradation, 








The lack of monomer and the presence of Siky and 


CO, are noted. A very small quantity (<1%) of 


other silicon derivatives were also observed but were 
not identified. These results appear to indicate a 
reaction of the type: 

CF,—CF, - +O, products. 

The primary products may consist of labile com- 
pounds, possibly CFO, that quickly decompose on 
the glass at the temperature of pyrolysis to give the 
observed compounds. There appears to be no attack 
on the polymer because this would be expected to 
result in a darkening of the polymer and or a drop 
in molecular weight to give a residue of high fluidity, 
at the temperature of pyrolysis. 

When present at high concentrations in the thermal 
decomposition of hydrocarbons, nitric oxide has been 
shown [5] to catalyze the rate of decomposition. At 
the high concentrations used here, the nitric oxide 
also seems to have an over-all cataly tic effect. The 
polymer darkened slightly initially, and as conver- 
sion increased, the residue became clear and very 
fluid. This, plus the fact that the rate starts to level 
off, would seem to indicate that not only an induced 
initiation by some transfer mechanism is” taking 
place, but also some product is being formed that 
stabilizes the degrading free radicals or polymer. 

In the experiment with hydrogen sulfide the poly- 
mer quickly blackened and became fluid at its tem- 
perature of pyrolysis. During the course of this 
experiment the apparatus became coated with sulfur. 

As the apparatus was etched from other experi- 
ments, it was not known definitely whether or not 
hvdrogen fluoride was produced here; it presumably 
was. These results would seem to indicate induced 
initiation and defluorination. 

With sulfur dioxide the polymer became black and 
gelatinous. Again, this is indicative of induced de- 
composition and defluorination with no inhibition, 
The increasing rate of degradation as shown in figure 
4 supports this. 


4.2. Inhibition Effects 


In figures 5, 6, and 7 are plotted the weight-loss- 
versus-time data for the thermal decomposition of 
polytetrafluoroethyvlene under hvdrogen, chlorine, 
and toluene. Carbon tetrachloride also showed 
good inhibition initially, but as at later stages, the 
inhibition was much less effective than with the 
above three compounds, the system was not studied 
in detail. 

It will be immediately noted on comparing these 
results with those obtained for the pyrolysis of poly- 
tetrafluoroethylene under nitrogen, figure 2, that the 
characteristics of the reaction have completely 
changed in that all the curves for the gas studies have 
a definite inflection point. Consequently, when the 
rate is plotted against conversion, figures 8, 9, and 
10, it is seen that all the graphs have definite maxima 
and that initially the rates are much lower than those 
for the vacuum work, figure 1. These curves are 
tvpical of those obtained for mechanisms with large 
random components [2]. Random initiation or very 
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frequent transfer with short propagation and termi: 


nation usually give rise to such results. 


It has beer 


recently demonstrated that the thermal decomposi- 
tion of polymethyvlene follows such a pattern [6 
This type of mechanism leads to a very fast drop i 
molecular weight, and this has been noted in these 
experiments, because the polymer quickly became 


very fluid. 


The experiments with chlorine as the inhibitor did 
not show any consistent color change in the polymer 
Also, an analysis of the polymer after SO-percent 
degradation showed a chlorine content of 0.6 percent 
by a microcombustion method. All this is compatr 
ble with the idea that the mechanism of the observed 
inhibition involves the combination of polymer radi- 


eals with the chlorine atoms produced, 
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of the gaseous products, table 1, is very complicated, 
consisting mostly of silicon, halogen, and carbon- 
containing compounds, ‘These again would seem to 


RATE OF VOLATILIZATION , PERCENTAGE OF ORIGINAL SAMPLE PER HOUR 
7] 








result from the decomposition of labile primary ' of qee"s 
products on the glass. Co onD 
In the experiments with hydrogen there is, besides a 6 
the molecular weight drop, a blackening of the 
polymer. At very high conversions (about 65%) IF 
the residue becomes hard. Concurrent with all this 
is an etching of the glass tube. An analysis of the 0 n 1 
gaseous products, table 1, shows that some monomer a: 3 2 2 oe 6S Um hUmhhUm 
is produced, as well as carbon dioxide and silicon VOLATILIZATION , PERCENT 


tetrafluoride. These results again indicate random 
degradation, but with the added feature of defluorina- 
tion by hydrogen. 


Figure 10. Rate of volatilization of polytetrafluoroethylene 
under toluene diluted with nitrogen as a function of conversion 
at several temperatures. 
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The experiments with toluene give results similar 
to those with hydrogen. However, defluorination is 
much more prevalent and may actually be rate de- 
termining, because the residue, after becoming liquid, 
soon hardens to a brittle mass at about 30-percent 
conversion. The analysis of gaseous products is 
given in table 1. Here again there is some monomer 
formed, as well as silicon tetrafluoride and carbon 
dioxide. 
may mean relatively less effective reaction of the 
vapor with decomposing polymer free radicals. 


5. Discussion 


In any kinetic discussion some idea of the various 
bond strengths is valuable. However, the carbon- 
carbon-bond strength in fluorocarbons is in doubt. 
A value of 124 keal for the carbon-carbon bond in 
hexafluoroethane has been reported by Dibeler, 
Reese, and Mohler [7] as obtained in electron-impact 
experiments. However, a recent, separate deter- 
mination of the ionization potential of the per- 
fluoromethyl radical by Lossing [8], using the 
electron-impact method, has lowered this value for 
hexafluoroethane to 95 keal. 

An analysis of the best thermochemical data on 
the subject would seem to favor a value about half 
way between these two extremes. Using the high 
and rather generally accepted value of the heat of 
sublimation of graphite, i. e., LZ (C)=170.6 keal, 
together with the heat of formation of carbon tetra- 
fluoride, 220 keal [9, 10], an average carbon-fluorine 
bond of 116 keal is obtained for this compound. 
Now in hydrocarbons the average carbon-hydrogen 
bond strength in ethane is about 8 percent less than 
those in methane. If in fluorocarbons the carbon- 
fluorine bond strengths drop off by approximately 
the same percentage in going from carbon tetra- 
fluoride to hexafluoroethane, then the average 
carbon-fluorine bond strength in hexafluoroethane 
should be about 107 keal. This, taken with the 
heat of formation of hexafluoroethane of 303 keal, 
as obtained by Duus [9], gives a carbon-carbon bond 
strength for this fluorocarbon in the neighborhood 
of 112 keal. 

If the interpretation of the kinetic data along the 
lines of the existing theories is correct and poly- 
tetrafluoroethylene does decompose by end initiation 
with the average carbon-carbon bond strength being 
close to that in hexafluoroethane, e. g., 110 keal, 
it should be possible to deduce the strength of the 
bond broken in the end initiation (see section 1), i. e., 
K,. Taking the only value available for the molecu- 
lar weight of polytetrafluoroethylene, i. e., 200,000 
[11], and realizing that for end initiation to have a 
predominant effect on the kinetics, this rate of initia- 
tion must be 10 times faster at about 500° C than 
that for random initiation, it is found that there 
must be a difference of at least 15 keal between the 
strength of the weak end bond and that of the 
average bond. This gives the value of 95 keal, 
which was used for end initiation. 

In the experiments described it appears that 
the various gases have in some cases altered the 


The somewhat larger amount of monomer 
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Figure ll. Activation energy slopes for the thermal degrada 
tion of polytetrafluoroethylene under hydrogen, chlorine, ang! 
toluene. ; 


mechanism of the degradation of polytetrafluoro. 
ethylene. The that act as inhibitors hay 
introduced a random quality to the over-all degrada 
tion. 

In this random type of degradation as ascer 
tained from the kinetic curves (fig. 8, 9, and 10 
the maximum rates are the most dependable and th 
most significant from which to calculate energies @ 
activation. The maximum rate of volatilizatio, 
has been related [2] to the specific rate constant for 
purely random reaction as follows: 


[dC Ja nash 


gases 





where (is conversion, ¢ is time, & is a compositi 
specific rate constant for the over-all reaction, | 
is the base for the natural logarithm, and Z is thy 
critical chain length for vaporization (i. e., th 
size for which rate of decomposition equals rate ¢ 
vaporization) and is here assumed to be a constant, 
In terms of elementary chain mechanism and _ th 
total radical concentration FP, 


k=k, +k R=k;(I/k,)'”? (when &, is small), 
where J is the rate of initiation of active centers, f 
' the rate constant for transfer, and 4, the rate con 
| stant for termination. 

Taking the maximum rate from the experimenta 
curves as proportional to the specific rate constan 
and making Arrhenius-tvype plots of the data 
straight lines are obtained, as shown in figure 1! 

From the slopes of these lines energies of active 
tion are obtained. In the case in which chlorine 
was the inhibiting gas, a value of 100 +4 keal® 
obtained. With hydrogen, a value of 118 +4 kea 
is obtained, and toluene gives 64 +5 keal. 

The experimental results and theory can ap 
parently be brought together if a mechanism of th 
following type is assumed, where, beside the thre 
elementary reactions mentioned above, there is 
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(4) Initiation. 


k 
R—C F,—R+A-—> AF+P+RCF;- (5) Transfer. 
ks 
R—CF,—CF,—R+A; — R—CF=CF—R+2AF _ (6) Transfer. 
ky” . 
R—CF,-+A-—+ R—CF,—A (7) 
ka?" on —_— 
R—C F.-+A;—- R—CF,—A+A (8) @ Termination. 
tee 
A-+A-+X — > A+ X (9) 


is the inhibiting gas, P and R—CF,- are 


lere Ay 
= polymer and polymer radicals, respectively, and 
X is the third body required for this reaction. 

In the experiments with chlorine the low bond 
energy of this molecule should favor the formation 


of these atoms as the potential chain carriers. As 
the concentration of these atoms would be relatively 
high, the possibility of a chain transfer by reaction 
5 and a molecular-weight drop would be increased, 
having an over-all effect of inducing the polymer 
decomposition. The lack of monomer in the gases 
evolved would seem to indicate that reaction 2 is 
suppressed, presumably owing to the termination of 
the polymer radicals by reactions 7 or 8, resulting in 
inhibition. The over-all chain termination would, 
then, seem to be by reaction 9, because the chlorine 
atoms would constitute the largest source of chain 
carriers. As the polymer is colorless throughout the 
degradation, reaction 6 would not likely be opera- 
tive, as that would be expected to cause darkening 
of the polymer. 

With this mechanism the over-all energy of ac- 
tivation would be expected to be related to those of 
the elementary steps as follows: 


E=}E,+E,—-}E" 
57 : 
100=—+E;—0 


k,=71.5 keal. 
Here 57 keal is the dissociation energy of the chlorine 
molecule [12]. This rather high value for 4; would 
seem to be the energy necessary for a chlorine atom 
to abstract a fluorine atom from the polymer, which 
reaction must be assumed in order to explain the 
high maximum rate and the molecular-weight drop. 
When hydrogen is the inhibitor, the situation 
is slightly different. It is thought that the chain 
initiation step is the breaking off of a volatile frag- 
ment from the ends of the polymer, as shown by 
reaction 1, followed by unzipping of the polymer 
radical to give some monomer, reaction 2 (see 
table 1), but finally terminating by reaction 8. 
This mode of initiation of the reaction is assumed 
in place of hydrogen dissociation because the bond 
strength of hydrogen is higher than the 95 keal of 
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| this reaction. Chain unzipping would be inhibited 
by reaction 8, and the hydrogen atom so formed 
probably induces further decomposition in the 
polymer by reaction 5. <A side reaction such as 
reaction 6 would account for the darkening of the 
polymer but, due to the small mass of material 
vaporizing in this way, would probably not drasti- 
cally affect the over-all rate. Final termination of 
the over-all chain reaction would seem to be most 
likely by reaction 7 if a steady state is truly obtained. 

In this inhibition under hydrogen, the over-all 
energy of activation for this random process is 
broken down as follows: 


E=}F,+F,—}E/ 
95 . 
113 _ 9 4 hk,—0 
k,=70.5 keal. 


These values obtained for the energy of activation 
for transfer, i. e., £3, are high and constant within 
the experimental error of the experiment, e. g., 
71 +4 keal. As this step probably comprises two 
steps, one the abstraction of fluorine by a radical, 
and the other the breaking up of the polymer radical 
in a depropagation step, this high and apparently 
constant value is probably due to the propagation 
step £). 

In the case of the toluene inhibition, an initial 
drop in molecular weight and the defluorination 
can be explained by the above reaction scheme, as 
in the case of hydrogen and chlorine. Here, forma- 
tion of hydrogen atoms and benzyl radicals may be 
the most important initiation step (reaction 4), 
which in time should lead to chain transfer, reaction 
5, and hence induced decomposition of the polymer. 
Termination by reactions 7 and 9, where benzyl 
radicals are important, would end the chain reaction. 
The observed inhibition could be explained by 
reactions 7 and 8. The presence of the defluorina- 
tion by reaction 6, as evidenced by early hardening 
of the polymer, complicates the over-all rate. 
Therefore, it is not possible to identify the energy 
of activation obtained from the toluene experiments, 
i. e., 64 +5 keal, as a single entity or a function of 
elementary energies of activation. 
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Abscissas and Weights for Gaussian Quadratures 


of High Order 


P. Davis and P. Rabinowitz 


Weights and abscissas are presented for the Gaussian quadrature rules of order n 
These constants were computed on Standards Automatic Computer by 


20, 24, 32, 40, 48. 


16, 


the method described and have passed a number of checks with about 20 places of decimals. 


Values of the weights and abscissas are also available for n 


1. Introduction 


Gaussian quadrature rules have hitherto been 
computed up to the case n 
curacy [1]. 
Section of the National Bureau of Standards, fre- 
quent use has been made of the rule n=16 and even 
of halving and quartering the interval for increased 
accuracy. For this reason, it is felt that the con- 
stants for rules of higher order will prove to be of 
use in working with electronic digital computers. 

For further evidence of the practical utility of 
high-order rules, the reader may consult Hartree [2], 
Henrici [3], and Reiz [4]. Exact values of these 
quantities are also interesting in view of certain 
unsettled theoretical conjectures that have been 
made about distribution of the weights and abscissas 


[7]. 
2. Method of Computation 


We deal with the Gaussian quadrature rule of 
order n on the interval [—1,1}]: 


*+1 n 
| Aa)dz => in f(ren). (1) 
—] k=1 


The rule (1) holds exactly whenever f is a polynomial 
of degree <2n—1. The abscissas v,,(4=1,2,..., n) 
are the n zeros of the Legendre polynomial of order 
n: P(4¢,)=0, whereas the weights are given by the 
expression 
Ak, +1 in l 
’ 
k, Pn + (ry n)P al 2) a) 


Ary 


where ppt,—h, nF . are the normalized Legendre 


polynomials. See, e. g., Szegd [5, p. 47]. Making 
use of the relationship 
(—2)Pi(r)=nlP,_(7)—2rP,(2)], (3) 


we are able to derive from (2) the following alternate 
expression for the weights a,,, which is useful for 
computation: 


2(1 — a7.) 
= », 60 
InP, (Len) ]* 


SEs 


rr a |" (4) 


' Figures in brackets indicate the literature references at the end of this paper. 





16 with 15 place ac- | 
In the work in the Numerical Analysis | 





64, 80, and 96. 


To obtain a first approximation to the zeros of the 
Legendre polynomials, we make use of the following 
inequality derived by Szegé [6]: Let 


Tin =COS Oy» (k oa (5) 
Then 
Je P Ca jr e 
[(m 1/2)? e/4] 2 "<4 1/2) (k=1,2 ) 
(6) 


where j, (k=1,2, .) are the successive zeros of 
the Bessel function Jy (4) and e=1— (2/r)?’. 

The first 150 values of j), may be found in [8]. A 
preliminary computation of (6) with n=16 showed 
that the value of 6,, is closer to the left-hand bound, 
and that five or six decimal places can be secured 
intially by employing 


j 


((n-4-1/2)?+¢/4)'? (7) 


2; =COs 


as a first approximation to z,,. The value 2? 
was successively improved by using the Newton 
formula 


> ) 
r i+1) — r i _f jilte ) (8) 
en Pua) 
The derivative in (8) was computed from (3), 


whereas the Legendre polynomials themselves were 
computed from the recursion 


nP,(2r)- 


(2n—1)2P,,_,(4)—(n—1)P,_2(2) } (9) 
P(r) 


i. P(r) z. 
After the first approximation, the successive approxi- 
mations were computed in double precision, and a 
shutoff value of e=2~" was employed in the itera- 
tion (8). 
Although the abscissas and the weights are 
symmetric about s=0, all were computed inde- 
pendently. The starting values z;,(') for two 
symmetric points were different (cf. (7) ). This 
served as one check of the accuracy of the compu- 
tation. Additional checks were provided by com- 
puting the six quantities given in eq (10), page 37. 
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Abscissas 


0.5773502691 


0.8611363115 
0.3399810435 


0.9602898564 
0.7966664774 
0.5255324099 
0.1834346424 


0.9894009349 
0.9445750230 
0.8656312023 
0.7554044083 
0.6178762444 
0.4580167776 
0.2816035507 
0.0950125098 


0.9931285991 
0.9639719272 
0.9122344282 
0.8391169718 
0.7463319064 
0.6360536807 
0.5108670019 
0.3737060887 
0.2277858511 
0.0765265211 


0.9951872199 
0.9747285559 
0.9382745520 
0,8864155270 
0.8200019859 
0.7401241915 
0.6480936519 
0.5454214713 
0.4337935076 
0.3150426796 
0,1911188674 
0,0640568928 





89625764509 


94052575224 
84856264803 


97536231684 
13626739592 
16328985818 
95649804939 


91649932596 
73232576078 
87831743880 
55003033895 
02643748447 
57227386342 
79258913230 
37637440185 


85094924786 
77913791268 
51325905868 
22218823395 
60150792614 
26515025453 
50827098004 
15419560673 
41645078080 
33497333755 


97021360180 
71309498198 
02732758524 
04401034213 
73902921954 
78554364244 
36975569252 
88839535658 
26045138487 
96163374387 
73616309159 
62605626085 


n= 20 


n=24 


Weights 


1.0000000000 


0.3478548451 
0.6521451548 


0.1012285362 
0.2223810344 
0.3137066458 
0.3626837833 


0.0271524594 
0.0622535239 
0.0951585116 
0.1246289712 
0.1495959888 
0.1691565193 
0.1826034150 
0.1894506104 


0.0176140071 
0.0406014298 
0.0626720483 
0.0832767415 
0.1019301198 
0.1181945319 
0.1316886384 
0.1420961093 
0.1491729864 
0.1527533871 


0.0123412297 
0.0285313886 
0.0442774388 
0.0592985849 
0.0733464814 
0.0861901615 
0.0976186521 
0.1074442701 
0.1155056680 
0.1216704729 
0.1258374563 
0.1279381953 


0000000000 


3745385737 
6254614262 


9037625915 
5337447054 
7788728733 
7836198296 


1175409485 
3864789286 
8249278481 
5553387205 
1657673208 
9500253818 
4492358886 
5506849628 


3915211831 
0038694133 
3410906357 
7670474872 
1724043503 
6151841731 
4917662689 
1838205132 
7260374678 
3072585069 


9998719954 
2893366318 
1741980616 
1543678074 
1108030573 
3195327591 
0411388827 
1596563478 
5372560135 
2780339120 
4682829612 
4675215697 
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Abscissas 


0.9972638618 
0.9856115115 
0.9647622555 
0.9349060759 
0.8963211557 
0.8493676137 
0.7944837959 
0.7321821187 
0.6630442669 
0.5877157572 
0.5068999089 
0.4213512761 
0.3318686022 
0.2392873622 
0.1444719615 
0.0483076656 


0.9982377097 
0.9907262386 
0.9772599499 
0.9579168192 
0.9328128082 
0.9020988069 
0.8659595032 
0.8246122308 
0.7783056514 
0.7273182551 
0.6719566846 
0.6125538896 
0.5494671250 
0.4830758016 
0.4137792043 
'0.3419940908 
0.2681521850 
0.1926975807 
0.1160840706 
0.0387724175 


0.9987710072 
0.9935301722 
0.9841245837 
0.9705915925 
0.9529877031 
0.9313866907 
0.9058791367 
0.8765720202 
0.8435882616 
0.8070662040 
0.7671590325 
0.7240341309 
0.6778723796 
0.6288673967 
0.5772247260 
0.5231609747 
0.4669029047 
0.4086864819 
0.3487558862 
0.2873624873 
0.2247637903 
0.1612223560 
0.0970046992 
0.0323801709 


49481563545 
45268335400 
87506430774 
37739689171 
66052123965 
32569970134 
67942406963 
40289680387 
30215200975 
40762329041 
32229390024 
30635345364 
82127649780 
52137074545 
82796493485 
87738316235 


10559200350 
99457006453 
83774262663 
13791655805 
78676533361 
68874296728 
12259503821 
33311663196 
26519387695 
89927103281 
14179548379 
67980237953 
95128202076 
86178712909 
71605001525 
25758473007 
07253681141 
01371099716 
75255208483 
06050821933 


52426118601 
66350757548 
22826857745 
46247250461 
60430860723 
06554333114 
15569672822 
74247885906 
24393530711 
29442627083 
15740339254 
23814654674 
32663905212 
76513623995 
83972703818 
22233033678 
50958404545 
90716729916 
92160738160 
55455576736 
94689061225 
68891718056 
09462698930 
62869362033 


n= 32 


n= 40 


n= 48 





Weights 


0.0070186100 
0.0162743947 
0.0253920653 
0.0342738629 
0.0428358980 
0.0509980592 
0.0586840934 
0.0658222227 
0.0723457941 
0.0781938957 
0.0833119242 
0.0876520930 
0.0911738786 
0.0938443990 
0.0956387200 
0.0965400885 


0.0045212770 
0.0104982845 
0.0164210583 
0.0222458491 
0.0279370069 
0.0334601952 
0.0387821679 
0.0438709081 
0.0486958076 
0.0532278469 
0.0574397690 
0.0613062424 
0.0648040134 
0.0679120458 
0.0706116473 
0.0728865823 
0.0747231690 
0.0761103619 
0.0770398181 
0.0775059479 


0.0031533460 
0.0073275539 
0.0114772345 
0.0155793157 
0.0196161604 
0.0235707608 
0.0274265097 
0.0311672278 
0.0347772225 
0.0382413510 
0.0415450829 
0.0446745608 
0.0476166584 
0.0503590355 
0.0528901894 
0.0551995036 
0.0572772921 
0.0591148396 
0.0607044391 
0,0620394231 
0.0631141922 
0.0639242385 
0.0644661644 
0.0647376968 
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o n n cs ; ‘ 
te Se > Meakra=0, > enn =2 3 
eI k=1 t=1 
- n . 
de: ! 0, > Oxnxt =2/5, r (10) 
k=1 k=1 
2 n(n—l) B, 
Lin ~ 3/6 ‘ummm? Me 
= 2(2n—1) A, , 





where 


A,z"+B,2""'+...,n even. 


P,{2) 
These checks were all met to within 2 units in the 
90th decimal place. The first four checks in (10) 
were carried out on SEAC at the time of the com- 
putation, and the last two were made directly from 
the final tabulation. 

In the tables only the abscissas lying between 0 
and 1 have been listed. 


The authors thank the Hand Computing Unit of 
the Bureau’s Computation Laboratory for its 
assistance in checking these tables. 
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Behavior of Cements and Related Materials Under 
Hydrostatic Pressures up to 10,000 Atmospheres 


C. E. Weir, C. M. Hunt, and R. L. Blaine 


Compression studies were made on hardened cement pastes at 21° C and at pressures 
as high as 10,000 atmospheres. The effects of composition, age, moisture content, and 
water-cement ratio were studied. Compressibility of dried portland cement pastes appeared 
to proceed through a maximum with increased cure, the existence or location of the maximum 
being dependent on the water-cement ratio. Compression of both portland and aluminous 
cements increased with increasing evaporable water content. Portland cements with 
evaporable water contents in the range 17.5 to 32.4 percent exhibited a transition between 
2,000 and 3,000 atmospheres. No similar effect occured in aluminous cement. Experiments 
indicated that the transition was attributable to hydrated lime formed in the portland 
cement. Average compressibilities of dried paste were 3.1 10-° and 2.2 10-* atm be- 
tween 1 and 5,000 atmospheres and between 5,000 and 10,000 atmospheres, respectively. 


Compression studies were 
dry and wet calcium hydroxide. 


1. Introduction 


In considering studies of the pore-size distri- 
bution in hardened cement pastes by the mer- 
cury porosimetry techniques [1, 2, 3]'data are required 
on the magnitude of the compression of the cement 
when subjec ‘ted to hydrostatic pressures. Normally 
the compression of the cement paste would be 
considered negligible, but a study by mercury 
porosimetry of the extremely small pores believed 
to exist in dried cement paste would require rather 
high pressures at which the change in volume of the 
paste might well be of significance. 

Search of the literature revealed that while the 
elastic properties of concrete have been studied 
extensively, the measurements were conducted 
almost exclusively at low pressures and only rarely 
on neat cement [4]. However, the elastic properties 
of the rocks and minerals that may be used as 
aggregates in concrete have been investigated widely 
even at extreme conditions of temperature and 
pressure, primarily because of geological interest. 
The latter data have been compiled by Birch [5]. 

Because concrete is subjected to moderately high 
compressive stresses in many commercial applica- 
tions it appears that there is little direct evidence 
aside from successful usage —that the aggregate and 
the neat cement are compatible as regards com- 
pressibility. In order to obtain direct evidence on 
this question, as well as to evaluate the corrections 
required in high-pressure mercury porosimetry 
studies, preliminary data were obtained on the 
compression of hardened cement pastes with avail- 
able high pressure equipment. These data im- 
mediately showed evidence of a transition at ele- 
vated pressures in wet cement and indicated the 
desirability of a detailed study of the compressibility 
of cement itself, irrespective of any immediate sig- 
nificance of the data in direct applications. 





Figures in brackets indicate the literature references at the end of this paper, 





also made with wet silica gel, hydrated tricalcium silicate, and 


This paper contains the results of studies of the 
compression of neat cement pastes at pressures as 
high as 10,000 atm. Studies were made on the 
effects of cement type, water-cement ratio, water 
content, and length of cure. The compression 
characteristics of several pure materials of interest in 
connection with cement were also studied. 


2. Experimental Method and Preparation of 
Specimens 


2.1. Compression Measurements 


The apparatus, experimental technique, and 
method of calculating the data have been described 
in detail [6]. An experiment consists in forcing a 
leakproof piston into the smooth bore of a high- 
pressure vessel containing the specimen immersed in 

light petroleum distillate (Varsol). Readings of 
the depth of penetration of the piston are made at 
each 1,000 atm by means of a sensitive dial gage 
fastened to the piston, pressure in the vessel being 
measured by means of a calibrated manganin resist- 
ance gage immersed in the distillate. The apparatus 
was calibrated with steel bars having volumes 
similar to those of the specimens. The compressi- 
bility of steel has been measured by Bridgman [7]. 
The compression is defined as —AV/Vo, where Vo is 
the volume of the specimen at 1 atm, and AV is the 
change in the volume reckoned from the fiducial 
pressure of 2,000 atm; compressions at lower pres- 
sures, therefore, appear as negative values. All 
measurements described here were made with de- 
creasing pressures, and frictional forces at the pack- 
ing on the piston prevented the attainment of 
reasonable accuracy at pressures much below 2,000 
atm. Values of compression reported at 1 atm are 
not experimental but are obtained by successive 
approximations. Values of Vo were obtained by 
weighing in air before compression measurements 
and by “hydrostatic weighing in distillate after the 
measurements. 
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2.2. Preparation and Treatment of Specimens 


Test specimens were made by casting cement paste 
in paraflin-coated paper molds designed to yield a 
cylindrical specimen % in. in diameter and 4 in. long. 
Specimens were removed from the molds after 24 
hr and cured at 21° C over water in closed glass 
containers. Cured specimens were dried in a 
vacuum at room temperatures with a dry-ice trap 
interposed between the specimen and the pump in 
an apparatus similar to that described by Copeland 
and Hayes [8]. During the initial 24 hr the drying 
was carried out at a pressure of approximately 1 to 
2 mm of Hg. Following the initial 24-hr drying 
period, further drying was conducted at approx- 
imately 2 X 10°* mm Hg. Cylindrical test speci- 
mens prepared in this manner were immersed in dry 
distillate until evolution of displaced air appeared to 
cease and were then subjected to a short evacuation 
at a pressure sufficiently low to cause the distillate to 
boil at room temperature. The resulting vacuum- 
impregnated specimens were used for test. 

Certain materials studied could not be cast con- 
veniently in cylindrical form. Such materials were 
packed into a thin-walled steel cylinder, which was 
open at one end, having dimensions ': in. in diameter 
and 4 in. long. The material was impregnated with 
distillate, as previously described, prior to test. 


2.3. Materials Studied 


The materials studied included portland and 
aluminous cements, silica gel, tricalcium silicate, and 
calcium hydroxide. 


a. Cements 


The portland and aluminous cements studied in 
greatest detail had the compositions shown in table 1. 
The portland cement was denoted as cement A. 
This was a type-I (Federal Specification SS-C-—192a) 
cement ground from clinker, the composition of 
which is given in the table. Limited studies were 
also conducted on other portland cements meeting 
the requirements of Federal Specification SS-C-192a 


for types I, II, III, IV, and V. 
b. Silica Gel 


The silica gel studied contained 205 percent of 
evaporable water, based on the dry weight of the 
sample. An additional 5.4 percent of water was 
removed by ignition at 1,000° C for 15 min. After 
an overnight evacuation at 100° to 110° C, the 
surface area of the material by low-temperature 
nitrogen adsorption was 749 m?/g. This is a rather 
large surface area even for silica gel. 


c. Tricalcium Silicate 


The tricalcium silicate was prepared by heating 
and reheating of an intimately mixed paste of calcium 
carbonate and silica gel. 

d. Calcium Hydroxide 


The calcium hydroxide used was a reagent-grade 
product obtained commercially. 
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TABLE 1. Compositions of cements studied 
[- = 
Cement A®* (portland cement, type I) Aluminous cement 
Per- , Per- 
i centage, Calculated composition, ide centage 
Oxide by | percentage by weight Oxide | by 
weight | weight 
CaO... 65.13 | 3CaO-SiO2, 56% CaO... 36.8 
SiOz. 22. 62 3CaO0-AlLOs, 3% SiOz. 8.6 
AlOs 4.41 4CaO-AhO;- Fe203, 16°.) AlhOs 42.7 
FeO; 5.33 | 2CaO-SiO2, 23% ..| FezO3- 4.3 
SO; 0. 10 = ; SO; 0.3 
MgO 1.71 MgO 8 
Na20-. 0. 04 FeO 6.2 
K,O . 28 TiO. 0.4 
Loss on igni- | .37 . Gain on ig- aa 
tion.* nition, 
Insoluble resi- 4 Insoluble res- 1.7 
due. idue, 


* This is an analysis of clinker that was interground with different amounts ¢ 
gypsum and blended to 1.8 percent of SO3. Ignition loss of blended cement wa 
1.69 percent. 


2.4. Determination of Volume and Cement Conten; 
of Specimens 


The specimens were weighed in air before test and 
weighed in distillate after completion of the com. 
pressibility measurements. The volumes of the 
specimens at 1 atm, Vo, were calculated from the 
relationship 
- wW—W, 

Vo= ‘) 

=F ( 
where w is the weight of the specimen in air, w, is the 
weight of the specimen in distillate, and d, is the 
density of the distillate. 

The specimens were finally ignited in platinum 
crucibles at 1,050° C for 15 min, after first burning 
off the distillate in a flame. The cement content of 
the specimens, (’, was calculated from the relationship 

CU=w,(1+.), (2 
where uw; is the ignited weight of the hydrated speci: 
men, and 7 is the fractional ignition loss of the original 
cement. 


3. Results and Discussion 


3.1. General Considerations 


Prior to a discussion of the experimental results 
it is advisable to consider details of the experiments 
that affect the experimental accuracy and the rigor 
with which the experimental data may be treated. 

In many measurements on porous materials the 
failure to remove all air from the voids may introduc 
serious errors. However, in these studies the spec: 
mens were impregnated with distillate at reduced 
pressures. Under such conditions, if the whole void 


| 
: 


r 








volume contains air trapped at the reduced pressure 


the volume of the air at 1,000 atm will be quite small; 


Of still greater importance is the fact that the change’ 
in volume of the trapped air between 10,000 and 1,00 
atm will be expected to produce a negligible error it 
these measurements. Two other factors tend t 
further decrease this small error. First, a larg 
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quantity of air was displaced by distillate, and 
secon |, trapped air has been found to dissolve rapidly 
in the distillate at high pressure. Experimental 
evidence will be given to confirm the belief that any 
errors 1D compressibility arising from trapped air 
are negligible. 

Penetration of distillate into the small voids in 


these experiments was accompanied by a rapid de- 
erezse in pressure with time. This effect occurred 
only when the specimen was initially subjected to 
high pressures and was absent on subsequent com- 
pressions. The rate of penetration was not measured 
but was observed to decrease with time and, although 
considerable time was allowed to attain equilibrium 
at the maximum pressure, it was most probable that 
slow penetration occurred during the initial measure- 
ment in each experiment. 

Penetration of distillate into the specimens pro- 
duced an apparent decrease in the quantity of liquid 
in the pressure vessel. The disappearance of liquid 
effectively increased the required displacement of the 
piston, so that in most instances the maximum 
stroke was exceeded before reaching 10,000 atm. 
Generally, more distillate was added to permit at- 
tainment of 10,000 atm, but in some instances 
measurements were made only to the maximum pres- 
sure attainable on the initial compression. In most 
of the latter instances the maximum was 9,000 atm, 
but in one instance it was only 8,000 atm. 

The experimental error in these data is subject to 
some uncertainty. In general, with the techniques 
used here it has been concluded that the data may 
he reproduced to within +0.0004 in —AV/Vo5 (9). 
Such precision is readily obtained when the correction 
terms [6] are kept reasonably small. In the present 
studies the correction terms were frequently un- 
desirably large and were in general not subject to 
control because of the penetration of distillate into 
the specimen and variations in the real volume of 
the specimens. The actual volume of the specimen 
studied was likewise subject to some uncertainty 
arising from occasional loss of small fragments in 
handling prior to weighing hydrostatically. 

Although this compilation of known sources of 
error is long, most of the errors involved are quite 
small and in many instances may mutually com- 
pensate. It is believed that each measured value of 
compression in these studies may be considered 
reliable to within +0.0008. The differences be- 
tween data obtained with duplicate specimens may 
be somewhat greater than this and represent differ- 
ences in the specimens themselves. Therefore, 
some of the conclusions based on small differences 
between different specimens are subject to some 
uncertainty. 


3.2. Compression of Dried Cement Pastes of Different 
Ages and Water-Cement Ratios 


Measurements were made on dry cement A before 


mixing it with water and after hydration for various 
periods of time. These specimens were prepared 


S691S1 ™) } 





4] 


with nominal water-cement ratios 0.3, 0.5, and 0.7 to 
determine whether the distillate penetrated all the 
voids or whether a portion of the change in volume 
was due to compression of voids. The data on these 
specimens, which were all dried in a vacuum for 
week, are given in tables 2 to 5. 

All compressions as tabulated are based on the 
reference pressure of 2,000 atm, but may be caleu- 
lated from 1 atm by adding to each value of the 
corresponding column the (positive) compression 
given at l atm. At the head of each column the 
nominal water-cement ratios are given. The actual 
water content was probably somewhat less than that 
indicated by the ratios given, which, however, serve 
to designate specimens differing significantly in 
total void content. Below each water-cement ratio 
are given in order the nonevaporable-water content 
based on the weight of the original cement, the 
density of the specimen, and the volume of the 
specimen at 1 atm. 


TABLE 2. Compression of hardened cement A paste cured for 1 


week 
Nominal water-cement 
ratio _. a 0.3 O.5 0.7 
Nonevaporable water 
content, ‘, of 
original cement 12.1 11.1 11.1 10.8 
Density at 21° C,g¢ em 2.679 2. 706 2.711 2.712 
Volume of test: speci- 
men, ¢m S. S61 5, 869 5. 702 4. 737 
Pressure Compression, ~AV 1 
atm 
10, 000 0.0193 0. 0193 0.0176 0.0181 
4, 000 OL7S O176 O155 O176 
s, 000 O61 0157 OSA 0162 
7. 000 0136 0127 Olly 0136 
6, 000 oll4 . 0106 0095S o112 
| 
5, 000 0090 . O80 OO74 . 0090 | 
4, 000 me . 0051 0048 0060 
3, 000 OO31 . 0020 0020 0030 
2, 000 0000 0000 . 0000 0000 
1,000 : —. 0016 ; 
1 ~_ (WM) —. 005 —. 005 | =, 006 


Compression of hardened cement 
month 


TABLE 3. A paste cured for 1 


Nominal water-cement 
ratio O38 0.5 0.7 


Nonevaporable Wate 
content, ©; of original 
cement 15.2 16.0 15.4 
Density at 21° C, g em 2. 559 2. 567 2. 584 
Volume of test speci- 
men, em? 8.951 6. 550 ». 211 
Pressure Compression, —AV/ Vo 
atm 
10, 000 0. 0255 0. 0233 0.0217 
4, 000 0211 0214 Oly 
&, 000 O195 OSS O175 
7, 000 0169 O1LS5S . 0148 
6, 000 0136 0129 . 0125 | 
| 
5, 000 O107 0097 . OORY 
4, 000 OO77 O04 . 0061 
3, 000 0036 0025 0050 
2, 000 000 . 0000 . 0000 
1 —. OOS —. 006 —. 006 








Compression of hardened cement A paste cure d for 6 


months 


TABLE 4 


Nominal water-cement 


ratio 0.3 
Nonevaporabl water 
content of original 
cement 16.7 20. 1 32.8 
Density at 21° C, gem 2. 557 2. 455 2. 333 
Volume of test speci- 
men, en S87 7. 111 Hi. 33S 
Pressure Compression, ~All t 
atm 
10, 000 0. 0205 0.0224 0. 0276 
4, 000 ODS6 o2l4 o26y 
S000 O16S ole O46 
7, 000 OSS O16S 0217 
6, 000 oLoy ol4l OIM4 
5, 000 O0OSS OLS OLSS 
1,000 OO5S OOSS5 oO1ly 
3, 000 ool OOS OOo 
2.000 0000 O008 0000 
l On OOS ol 
TasLe 5. Compression of unhydrated cement A powder 


Density at 21°C, 


gem OS] * 
Volume of test 
specimen, em 7.687 
Compression 
Pressure , 
Al 
atm 
ho, ooo ool 
000 UD 
s. 000 oO1nso 
7. 000 O11S 
(HM) isu 
5 000) oy 
LLL (mes 
4. 000 moss 
2.000 (a 
1000 Oos2 
l On) 
* Density by ASTM method CISS-44 is 3.18 ¢ em 


In table 2, measurements are given for two speci- 
mens of water-cement ratio 0.5. The differences 
noted in data on these duplicates, particularly 
at the higher pressures, are believed to be larger 
than the experimental errors involved and, there- 
fore, are indicative of differences in the specimens 
themselves. , 

The data in table 2 indicate very little difference 
in the compressions of cement pastes of different 
water-cement ratio after 1 week of curing. This 
strengthens the argument that the values obtained 
represent compression of the hydrated cement 
and do not reflect the compression of voids, except 
possibly submicroscopic pores believed to be present 
in all pastes regardless of the water-cement ratios 
The differences in porosity of the specimens were 
quite evident from their behavior on initial com- 
pression. Time effects attributed to penetration 
of distillate into voids were observed, which were 
roughly proportional to the water-cement ratio, 
being extremely large in the specimens with a ratio 
of 0.7 and practically nonexistent) in specimens 
with ratio 0.3. 
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cement ratio on com pression. 


, Water-cement ratio 0.3; >, water-cement ratio 0.5; @, water-cement ratio 07 


On comparing the compressions for the differen 
curing periods shown in tables 2, 3, and 4, interesting 
differences may be noted. The data indicat 
that the compression initially increases with in 
creased curing and then may reach a maximum 
and decrease with further curing, with the time 
of cure required to reach a maximum being de 
pendent on the water-cement ratio. This behavior 
is illustrated in figure 1, where the compressior 
at 10,000 atm is plotted against the nonevaporable| 
water content. It will be noted that a well-define/ 
maximum shown for the lowest water-comealy 
ratio; a similar maximum is absent in the highest 
water-cement ratio, Whereas the intermediate water} 
cement ratio shows an intermediate behavior. — Ii 
must be noted that the effects of curing perio 
and water content shown in figure 1 are 2 
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men variability. 
however, the results suggest that 
laid down in the first few weeks of hydration é 
being reinforced during subsequent hydration. This! 
process is more marked at lower water-cement! 
ratios, because there is not only less water bu 
there is also less space for deposition of new gel {11 
At low water-cement ratio the rate of hydration & 
also slower in the later stages than in specimens ¢ 
higher water-cement ratio. Possibly the max: 
mum might be observed at all water-cement ratio? 
if hydration were continued long enough. 
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Figure 2. Effect of nonevaporable-water content on specific 


volume of hydrated portland and aluminous cements. 


, Portland cement; @, aluminous cement 


The compression of the parent unhydrated cement 
is given in the data of table 5. This material was 
investigated as a powder and was studied in the 


steel evlinder as deseribed previously. Although 
this material was investigated as a powder, the 


data are considered comparable with data obtained 
with evlinders. Evidence for the validity 
of this conclusion has been obtained in experiments 
comparing compressions of powdered and molded 
specimens of organic polymers [10]. 

Useful by -produe ts of the compressibility measure- 
ments are density determinations on specimens that 
have been evacuated and then impregnated with 
distillate at pressure of 10,000 atm. The values 
obtained should represent a rather uniform 
penetration of the pore structure by the distillate. 
Powers and Brownyard [11] have expressed the 
volume of the solid phase in hydrated cement 
as the volumes of the original cement and the non- 
evaporable water according to the relationship 


cast 


so 


3) 


V, C0. Wala, ( 
Where VV. is the volume of 
weight of the original cement 
of the original cement, w, the weight of nonevapor- 
able water, and ¢, the hypothetical specifie volume 
of the nonevaporable water. In terms of specific 
volume of the solid phase, an equivalent relation- 


the solid phase, ¢ the 
, ¢ the specific volume 


ship may be written as 

/ A Yal'ns (4) 
Where, in addition to terms previously defined, 
ry, is the specific volume of the solid phase, and 
rand y, are the weight fractions of cement and 


nonevaporable water, respectively. With the nec- 
essary relationship 
r.=1—y, (5 
eq (4) may be transformed into 
r Vet (Un—Ve)Yn (0) 
From this equation a plot of y, against r, should 


vield a straight line of slope (v,—r,) and intercept ., 
on the v, axis. Values of y, the fraction of non- 
evaporable water in the dried specimen (expressed as 
percent), and v, are plotted in figure 2. A linear 
relationship is obtained which by extrapolation gives 
the specific volume of unhydrated cement at zero 
water content. The experimental value for unhy- 
drated cement in figure 2 was obtained with a Le 
Chatelier flask (ASTM standard method C188—44), 
because the valer obtained after compression 
appeared to be high and an error in weighing was 
suspected. The nonevaporable water content of the 
portland cement specimen with the highest specific 
volume in figure 2 is also suspected to be in error, 
because it is greater than values for specimens pre- 
pared from the same cement and cured a year. Such 
an error might easily occur through the loss of a frag- 
ment of specimen which would be caleulated as a loss 
of water. The nonevaporable water content of this 
specimen is given in table 4 as 32.8 percent of the 
weight of the original cement and is also plotted in 
figure 1. Powers and Brownyvard give 20 to 25 per- 
cent, based on the weight of the cement, as a reason- 
able estimate for the nonevaporable water content of 
a specimen close to ultimate hvdration {13}. The 
linearity observed below about 20-pereent non- 
evaporable water in figure 2 offers strong evidence 
that the distillate penetrated the same domains of the 
structure irrespective of age or water-cement ratio. 
The exact dimensions of the finest pores in cement 
gel are somewhat conjectural but may be of the order 
of a few angstrom units in diameter [13,14]. It 
doubtful if distillate penetrated all these ultrafine 
pores. However, penetration undoubtedly occurred 
in all pores whose size or number could be varied by 
changing the water-cement ratio. 

From the slope of the line of figure 2 and the value 
of r., a value for the hypothetical specific volume of 
the nonevaporable water, 7,, was calculated. The 
value obtained was 0.87 ml/g. This value is higher 
than the 0.82 ml/g reported by Powers and Brown- 
vard from measurements of density of a number of 
cements by helium displacement {11] or the value of 
0.74 ml/g obtained by Copeland {12], using water as a 
displacement medium. This indicates that distillate 
did not penetrate the structure to the same extent as 
did helium or water. Another factor contributing 
to this difference may be the difference in the method 
of drving. These specimens were dried as ‘s- by 
4-in. evlinders, whereas the specimens of Powers and 
Brownyard [11,15] and Copeland [8,12] were dried in 
granular form. Somewhat more water may be 
moved in the latter form. 
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3.3. Effect of Evaporable Water on Compression of 
Hydrated Cement ° 


The effect of evaporable-water content on the 
compression of cement of a given age and water- 
cement ratio was most interesting. In these exper- 
iments six specimens of cement A mixed with a 
water-cement ratio of 0.5 were cured for | week. 
Following the curing period the specimens were dried 
in a vacuum for periods varying from 15 min to 1 
week before the testing. The specimen dried for 1 
week retained 10.9 percent of water, based on the 
weight of the cement used in preparing the specimen. 
Although more rigorous drying procedures might 
have removed a portion of this water, the 10.9 per- 
cent of water was arbitrarily taken to represent the 
ionevaporable-water content of all the specimens. 
The compression data obtained on these specimens 
are given in table 6, and the compression-pressure 
behavior is shown graphically in figure 3. In figure 
3 the compressions of the specimens with evaporable- 
water contents of 2.2 and 5.2 percent were only 
slightly greater than the driest specimen. The 
specimens with evaporable-water contents of 17.5 
percent or more, however, were highly compressible 
and the compression increased markedly with in- 
creasing water content. Although the increased 
compression extends over the whole pressure range, 
it will be noted in the three specimens of high mois- 
ture content that the curves are nearly parallel 
above about 5,000 atm. 


Effect of « va porable water on compression of hardened 
cement A paste (1-wee k cure 


TABLE 6 


Evaporable-water 
content ot orig- 
inal cement $2.4 ny 17 2 2.2 “ 
Density, gem 1. 951 2 100 2. 270 ) So 2 HO 2.735 
Volume of test speci- 
nen, em 10. SOS sey 7.874 H.3l OF | 7. 708 
Pressure Compression vit 
wm 
10. 000 0.0777 0.0219 0.0225 0.0182 
4000 0733 «0.06089 «(0.04605 O1US 204 ol72 
000 ange Ons 0426 O178 O1S2 Os! 
7.000 iW O579 O375 oO14s O1ss O126 
6,000 O5S9 O524 0323 O22 O127 ole2 
5000 O52s i442 272 Ongs oloo wo74 
ion 445 O376 216 oon O08 OO4S 
3, 000 OS26 0253 O130 oo2s ony OO1T 
2, 000 0000 0000 0000 0000 0000 0000 
1,000 O16S 0120 O127 
l 03 v2 “2 OOn oon on 
. eg ene . 
Compressibility—the slope of the curves—is less 


affected by moisture content above 5,000 atm even 
at high evaporable-water content. Below about 
4,000 atm a rapid change in compression is observed 
at the higher moisture contents, and the compression 
in this region of pressures is a strong function of the 
water content. Below 4,000 atm the sigmoid shape 
of the compression curve clearly indicates that a 
transition is occurring. The transition is not sharp 
and appears to occur principally between 2,000 and | 
3,000 atm. Large thermal effects are associated 
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with this transition and these are readily observed 
experimentally [9] in their effect on the pressure 
From the experimental observations of the thermal 
effects it concluded that the transition not 
occurring below 2,000 atm but that it persists to 
pressures as high as 4,000 atm. No experimental 
evidence was found for any comparable discontinuity 
at the lower moisture contents, and it is concluded 
that the phenomenon requires moisture. 

The transition is completely analogous to a phase 
transition produced isobarically by changing the 
temperature, i. e., freezing point. The fact that this 
transition occurs over a rather wide pressure interval 
is understandable because in a complex material 
such as cement many other substances may enter 
into the phase change. In this respect it is analogous 
to the freezing of an impure material. 

A different picture of the process may be obtained 
by considering the compression isobarically — in 
relation to the evaporable water content. Isobar 
at a few representative pressures are given in figure 
4. The broken lines in figure 4 represent compres 
sions calculated for ideal mixtures of water and dry 
paste and converge on the compression axis with 
the corresponding experimental isobars. The calev- 
lated compressions were based on the known com- 
pression of water [16] and dry paste (the specimen 
containing no evaporable water in table 6) and their 
respective volume fractions. Significant divergence 
of the caleulated and experimental curves appear 
at pressures of 3,000 atm and higher and approx:- 
mately 10 percent or more evaporable moisture 
At 1,000 atm the calculated and measured data are 
in substantial agreement over the whole moistun 
range. 
calculated and measured values increases rapidly 
with increasing moisture content but is relatively 
insensitive to changing pressure at a given moisture 
content at the higher pressures. These observatiors 
are consistent with the behavior to be expected from 
a transition centered between 2,000 and 3,000 atm 
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Figure 4. Calculated and measured isobars for portland 


cement as functions of evaporable-water content. 


Digits at right of each isobar denote the presence in kiloatmospheres 
, Measured; . calculated 


Water-vapor sorption measurements have been 
made on separate specimens of hardened paste pre- 
pared from the same cement and cured 1 week over 
water. If the adsorption of water vapor by hy- 
drated cement is a multilayer adsorption process, in 
the range of relative humidity over which = the 
Brunauer-Emmett-Teller (BET) equation applies, 
V, is the amount of water theoretically required 
tocover the total surface of the hydrated cement with 
a laver of water 1 molecule thick. A V,, value 
of 2.8 percent of evaporable water was calculated 
from such data by the BET equation [17]. From 
figure 4 it is noted that the correlation between the 
experimental and calculated compressions deviates 
appreciably above an evaporable water content 2 
or 3 times V,,. 

The data at low moisture contents are of interest 
in connection with questions concerning the nature 
of this water. Within the experimental errors no 
essential differences exist between calculated and 
measured compressions at low moisture contents as 
shown in figure 4. Although the experimental points 
at 5.2 percent of evaporable water are indicated as 
lying somewhat below the calculated curves, there 
issome doubt as to whether this is of real significance. 
However, all data shown in figure 4 are based on 
evaporable-water content, and all specimens are con- 
sidered to contain 10.9 percent of nonevaporable 
water. The compression curve for the specimen 
containing 10.9 percent water has been compared 
with a similar curve for unhydrated cement and also 
with a curve calculated for a mixture of unhydrated 
cement and 10.9 percent of water, assuming the 
specific volume of nonevaporable water to be 0.87 


ml/g. These results are shown in figure 5. Similar 
curves could be made from the data in tables 
2to 5. Although all the compression curves would 
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Figure 5. Calculated and measured compressions for portland 
cement as affected by nonevaporable-water content. 


@, Calculated (unhydrated cement and mixture of 10.9 percent of liquid water); 
, measured (hydrated cement with 10.9 percent of nonevaporable water); 
, measured (unhydrated cement 


not coincide as closely with the unhydrated cement 
as the specimen shown in figure 5, all results would 
show that the experimental compressions were much 
less than those calculated for corresponding mixtures 
of unhydrated cement and water. Thus compressi- 
bility measurements confirm what would seem evi- 
dent from other considerations, namely, that non- 
evaporable water is accommodated differently in the 
structure of hydrated cement than evaporable water. 

Before describing the experiments that were con- 
ducted to ascertain the source of the transition ob- 
served in hydrated portland cement containing 
large amounts of evaporable water, it is of interest to 
consider certain possible explanations for the phe. 
nomenon. Water is known to freeze to ice VI at 
approximately 8,700 atm at room temperature [18]. 
Liquid water in cement will exist in a solution which 
will be expected to freeze at still higher pressures. 
A marked reduction of this freezing point by some 
6,000 atm is not expected. It has been proposed, 
however, that adsorbed moisture and water present 
in fine capillaries are subject to extreme surface 
forces and, therefore, normally exist under high 
hydrostatic pressure [19]. Such concepts have fre- 
quently been advanced in connection with theories of 
bound water which will not be reviewed here. A 
small increase in external hydrostatic pressure on 
such tightly bound water might exceed the freezing 
pressure and produce the transition observed here. 
However, such an explanation should lead to obser- 
vations of a transition even at low moisture contents 
at which it appears to be completely absent. In 
addition, no similar effects were observed in studies 
on fibrous polymers containing as much as 30 percent 
of moisture [20]. It has been suggested, also, that 
the high hydrostatic pressure might rupture small 
particles of cement as a result of nonhydrostatic 
forces at grain boundaries, thereby producing an 
apparent transition. However, all data reported 
here were obtained on decreasing pressure and in 
these studies measurements may be made only 








in this direction because of penetration of distillate 
into the voids upon compression. However, additional 
observations indicated that the transition was com- 
pletely reversible, and it is difficult to conceive of a 
reversible but destructive physical phenomenon. 
Furthermore. as later data will show, a cement with 
no detectable transition has been studied. This 
fact alone essentially eliminates the forementioned 
explanations for the transition. 

From the data already presented it appears that 
some hydrated material in cement produced the ob- 
served transition. The wide pressure range of the 
transition is attributed to the fact that hydrated 
cement is not a pure substance 


3.4. Source of the Transition in Hydrated Cement 
Containing Large Amounts of Evaporable Water 


A number of experiments were performed to study 
the transition in detail and to ascertain the material 
responsible for the phenomenon. In order to de- 
termine whether the transition disappeared per- 
manently once the sample was dried, or whether it 
could be regenerated by rewetting, the following ex- 
periment was performed with three specimens pre- 
pared from cement A, nominal water-cement ratio 
0.5, cured 1 week. The first specimen was tested 
immediately following the l-week cure, with most 
of its evaporable water intact. The other two 
specimens were dried in vacuum for 1 week and one 
was tested in the dried condition. The last specimen 
was rewet in liquid water and tested immediately. 
The results obtained in these experiments are shown 
graphically in figure 6. It is apparent that the 
transition appearing in the original specimen is 
eliminated by drving and regenerated on rewetting 
the dried specimen. Whatever the source of the 
transition, it must arise from a material formed 
rapidly in the presence of water. 

Compressibility measurements were made on 
hardened aluminous-cement pastes containing differ- 
ent amounts of evaporable water. This was done 
with the idea that some structural similarity might 
exist between the hydrated tricalcium aluminate of 
portland cement and some of the hydrated com- 
pounds in aluminous cement paste [21]. No transi- 
tion was observed with hydrated aluminous cement 
containing large amounts of evaporable water. 
These data will be discussed in the next section. 

Compression measurements made on silica gel 
containing a large amount of evaporable water 
showed no transition within the pressure range in- 
vestigated. A specially prepared cement containing 
no gypsum was hydrated and observed to exhibit a 
smell transition on compression. These data, to- 
gether with those on the silica gel, are only qualita- 
tive and are not given in tabular form. The graphi- 
cul data, however, are shown in figure 7. It is con- 
cluded that neither gypsum nor calcium sulfoalu- 
minates were involved in the transition. 

Further studies were made using a specimen of 
tricalcium = silicate mixed with approximately an 
equal weight of water and cured for 1 week. Analy- 
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sis of part of the specimen showed that it contained 
9.9 percent of nonevaporable water, expressed ® 
percentage of the tricalcium silicate. The cOomMpressiol 
data showed a very strong transition between 2,00 
and 3,000 atm, as shown in figure 7. The transition 
here is much sharper than observed in the cements} 
but occurs in the same pressure interval. On drving 
and remeasuring the test specimen, the wransitio} 
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As free Ca(OH), is one of the hydration produets 
of tricalcium silicate and other cement compounds 
22], a final experiment was performed on Ca(OH), in 
the wet and dry states. A very strong transition was 
observed in the wet lime, which was completely 
absent in the dried material. The graphical data on 
wet Ca(OH), are given in figure 7. It will be noted 
that the transition in this material occurs at ap- 
proximately the same pressure and is noticeably 
sharper than that in the tricalcium silicate. 

Although the transition in wet cement is extended 
over a pressure interval, the major portion of the 
volume change occurs at essentially the same pressure 
at which the transition takes place in wet Ca(OH),. 
The chance of finding two different materials having 
similar transition temperatures at 1 atm is quite 
small. The chance of finding different materials 
having the same transition pressure at a given tem- 
perature is remote because pressure-induced transi- 
tions are less common than those of thermal origin. 
Itmay be concluded, therefore, that the discontinuity 
in volume observed in hydrated cement contain- 
ing appreciable evaporable water erises from wet 
Ca(OH)., which ts liberated in the hydration of port- 
land cement. The svstem Ca(QH),.-H,O has been 
investigated under pressure in additional experiments 
reported elsewhere, and it has been concluded that 
the transition arises from the formation of a hvdrate 
of Ca(OH), [23]. 


3.5. Aluminous Cement 


The aluminous cements occur in a different region 
of the CaQ-SiO,-Al,O, svstem from portland cements 
22). The compounds 5CaQO-3ALO,;, CaQ-3Al,Os, 
and 3CaQ-5ALO, are believed to constitute the 
major portion of aluminous cements instead of the 
compounds 2CaO-SiO,, 3CaO-SiO,, 3CaQO-ALO,, and 
4CaO-ALO,-Fe.O,, which are believed to exist’ in 
portland cements [22, 24]. Aluminous cements 
hydrate more rapidly than portland cements, with a 
corresponding more rapid evolution of heat. Free 


lime is not one of the final end products of hydration 
in the aluminous cements, as in portland cement 
[25] 
at. 


Studies made on aluminous cement included the 


effects of evaporable water and length of cure. 
Compression data on these specimens are given in 


tables 7 and 8 and are shown in part in figure 8. 
The specimen of highest evaporable water content in 
table 7 received no drying treatment; the driest 
specimen was dried in vacuum for 1 week, and had a 
nonevaporable-water content of 29.2 percent. In 
figure 8 it is observed that the compression increases 
with evaporable-water content, but from observa- 
tions during the measurements it was evident that 
there was no transition similar to that observed in 
the portland cements. In general, it was found 
that the compressibility of the hydrated aluminous 
cement was somewhat higher than that of the 
hvdrated portland cement, but part of this differ- 
ence may arise from a somewhat higher nonevapor- 
able-water content in the aluminous-cement pastes. 


TABLE 7. 
hardened aluminous-cement paste cured for 1 weel, 


Eff ct of ‘ va porable -water content on compression of 


E vaporable-water con 


tent, ©, of original 
cement 21.2 Is. 1 15.8 oo 
Density, g/em 2. O92 2.131 2.173 2. 390 
Volume of test speci 
men, em YW. SU4 ¥. t42 ¥. ISO 7.327 
Press! re Compression, —~AWV 
atm 
10, 000 0. 0436 0. 0373 0. 0319 0, 0202 
4.000 402 O340 O20 ODSs6 
s, 000 O34 O307 (257 O167 
7, 000 0321 0267 0220 O42 
6, 000 0267 0218 O1S2 O116 
5, 000 0211 O16S o142 oul 
4,000 O14s O1ls O1lde2 OO8L 
3, 000 O06 oos2 OST 0032 
2, 000 0000 O000 0000 Q000 
1, 220 OO6S A —_ (004 : 
| ol4 O12 —.01 —, 006 
TABLE 8. Compression of hardened aluminous-cement paste 
Cure - 1 «iuy l month 74e month 
Nonevapol ible - Water 
content, ©; of original 
cement. 23.7 33.2 37.5 
Density at 21 aa 
g/cm 2. 446 2.318 2. 263 
Volume of test speci 
men, ¢m 7.027 7. 708 8. 348 
Pressure Compression, —~ AV) Vj 
atm 
10, 000 0.0136 0). 0385 0. 0208 
4, O00 0122 0350 O24 
8.000 O12 0317 0254 
7, 000 OOs4 O2S1 O22) 
#8, O00 Ons o246 O1S2 
5.000 (M47 O20 O13S 
on) OOS Olas OOu4 
3. 000 OOLS OOSS O42 
2 O00 ooo ooo ooo 
1, 000 Wy one 
l OOS O15 O00 








The measured specific volumes of the specimens 
of aluminous-cement pastes have been plotted in 
accordance with eq (6). The value for the unhy- 
drated cement was obtained with the Le Chatelier 
flask because no compression measurements were 
made on this material. The graph of these data is 
shown in figure 2 and again indicates a linear rela- 
tionship. The slope of the line, which is somewhat 
less than that found for portland cement, vields a 
value of 0.75 ml/g for the hypothetical specific 
volume of the nonevaporable water as compared to 
the 0.87 ml/g found for portland cement by the 
same method. This difference that on 
hydration the system, aluminous cement plus water, 
undergoes a greater contraction than portland 
cement. Lea and Desch [22] have cited the work of 
Le Chatelier [26] and Gessner [27] in support of a 
similar conclusion. 

A graph of the isobars of aluminous cement 
similar to figure 4 for the portland cement is not 
included because the number of specimens was 
limited. However, experimental compression values 
were lower than values computed on the assumption 
that the compressions of dried paste and evaporable 
water were additive. This is quite different from 
moist portland-cement paste and adds support to 
the view that hydrated aluminous cement under- 
goes no transition like that produced in hydrated 
portland cement, or if any such transition occurs 
while the paste is under pressure, it is greatly 
reduced in extent. 

According to the data in table 8 the compression 
increases with hydration up to 1 month, but it is 
less than the maximum at 74: months. In other 
words, aluminous cement also shows a maximum 
similar to the ones shown in figure 1. The initial 
increase in compression is considerably greater than 
was found in the portland cement, but at least a 
portion of the effect may arise from the greater 
nonevaporable-water content in aluminous cement. 


suggests 


3.6. Compression of Different Types of Portland 
Cement 


In these studies measurements were made on 
hvdrated cements of different types mixed with a 
water-cement ratio of 0.5. Results obtained with 
specimens cured 1 week and vacuum dried 1 week 
are given in table 9. Measurements on specimens 
cured 1 month and longer are given in tables 10 and 
11. These cements showed differences in rate of 
hydration, as evidenced by their different moisture 
contents after drying. However, the correlation 
between degree of hydration and compression was 
only fair, which suggests other factors also play a 
part in determining the behavior of hydrated cement 
under pressure. Hydrated type I] cement was con- 
sistently a little less compressible than the others. 
the types Il and V 
cements were measured at three ages so that plots 
similar to the ones shown in figure 1 could be made. 
Such plots are not given for these specimens, but it 


The compressions of only 





| is evident that the type V cement had a maximum 
similar to those shown in figure 1, and the type IT 
merely showed curvature. 


TABLE 9. Compression of hardened pastes of various types of 
cements (l-week cure) 
Type of cement I II Ill IV \ 
Nonevaporable wa- 
ter content, “> of 
original cement 15.0 11.9 17.5 12.0 10.8 
Density at 21 J 
gem 2. 597 2. 683 2. 524 2. 697 2. 7 
Volume of test speci- 
men, ¢m 6. 422 5. 745 6. 538 6. 271 Wy 
Pressure Compression, ~All 1 
atm 
10, 000 0. OLS9 0.0211 O.O170 0.0170 
9. 000 O77 0.0121 o1g2 O160 Ose 
s, O00 OSS O07 O167 O1sY9 O1sy 
7, 000 O1S6 Oo! ol44 o1ly oll 
6, 000 oll O77 OLIS OO95 aogy 
5, 000 OOS6 Oos52 OOSS wo72 ord 
4, 000 Oo54 OOS5 OosY Lise} woe 
3, 000 0033 oo12 OOS0 OOS! oor4 
2, 000 0000 0000 OOOO Oooo (wh 
1, 200 - OO1S oOos2 ool OOR0 
l OO) -. O05 On) ws OOS 
TABLE 10. Compression of hardened pastes of various types of 
cements (l-month cure 
Type of cement II I\ \ 
Nonevaporable water 
content, ©; of original 
cement 16.9 Is. 2 128 
Density at 21° C, g/em 2.518 2. 525 2. 657 
Volume of test speci- 
men, em 6. 769 6. 824 ‘ ‘ 
Pressure Compression, AYv/' 
arm 
10, O00 0. 0207 OO, 
4, 000 ord Olds 
Ss, 000 O.O1sl O17 ol7o 
7, 000 O11 olay oOl45 
6, 000 O105 O124 OLIS 
5. 000 oovl mousy ol 
4,000 (O57 008g OosY 
3, 000 my OOSh oy 
2. 000 O000 OOOu OOO) 
l OO oo; on 
TABLE 11. Compression of hardened pastes of various types of | 


1 


cements (7'o-month cure 


Type of cement I I! \ 
Nonev por ible water 
content, ©; of original 
cement 25. 9 21.3 20 
Density at 21° C, g/cm 2. 367 2.4 2. 481 
Volume of test speci 
men, em 7. 687 7. d00 6.782 
Pressure Compression,—A) 1! 
atm 
10, 000 0. 0237 0.0177 0.0190 
4, 000 WZ OL7S Osh 
8, 000 0205 O61 olen 
7,000 OL76 ol44 ore 
6, 000 ~0146 O120 old 
5, 000 .O114 O0S8U 0006 
4,000 OOD) . OORT ou 
3, 000 0046 oo On38S 
2.000 O000 0000 . 0000 
1, 000 0045 0089 OO1T 
1 OOS 00H OO) 
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3.7. Compressibility and Elastic Constants of Cement 


The coefficient of compressibility, 8, is defined by 
meals Of the relationship 
8 dV/VidP, (7) 
where Vy is the volume at 1 atm and dV/dP the in- 
finitesimal change in volume divided by the infini- 
tesimal change in pressure. Values of 8 can be 
obtained as the slopes of the —AV/V, versus P 
eurves shown. For most of the dried cements 8 
decreases with increasing pressure only slightly and, 
therefore, average compressibilities above and below 
5,000 atm have been calculated from the equation 
8 AV/V,A4P, (8) 
where AV’ represents the change in volume between 
the maximum or minimum pressure and 5,000 atm, 
12 


TABLI -lverage compressibility of hydrated cements after 


vacuum drying for 1 week 


Average compressibility, l0-%atm 


Curing period 


lto 5.000 atm 45,000 to 10,000 atm 


{ 


Cement A: Water-cement ratio 0.3 
l-week 3.0 21 
l-month 3.7 3.0 
t-montl 2.9 2.6 
Cement A: Water-cement ratio 0.5 
l-week 2 6 92 
l-month $1 2.7 
temontl 4.4 te 
Cement A: Water-cement ratio 0.7 
l-week 30 18 
l-month 3.0 26 
temonth 5.2 24 
Pype I cement: Water-cement ratio 0 
l-week 2.9 2. 1 
7) »month 4.4 2.5 
ype Il cement: Water-cement ratio 0.5 
]-week 20 Ls 
l-mmonth 3.0 20 
7'.-mont! 4.0 1S 
Pyvpe IIL cement: Water-cement ratio 0.5 
l-week 3.0 20 
| vy" IV cement: Water-cement ratio 0.5 
l-week 24 20 
l-month 3.1 , 
Type V cement: Water-cement ratio 0.5 
l-week 2. 5 LY 
l-month $0 28 
7) month $1 LY 
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TABLE 14 


and AP the corresponding change in pressure. 
Average compressibilities calculated from eq (8) 


are given for a number of cements in tables 12 to 14. 

Table 12 contains data for dried specimens and 
shows the decrease in average compressibility with 
increasing pressure. Average compressibilities be- 
low 5,000 atm show a maximum with increasing 
curing time only for hydrated cement A, of water- 
cement ratio 0.3, whereas the corresponding data 
above 5,000 atm exhibit the maximum for all water- 
cement ratios for which the comparison could be 
made. Table 13 contains corresponding data for 
specimens of various evaporable-water contents. 
The marked difference in compressibilities noted 
above and below 5,000 atm at the higher moisture 
contents is directly attributable to the transition 
occurring in these specimens. Where such a dis- 
parity in compressibility is effected by pressure 
changes of this order of magnitude, it is most prob- 
able that a transition of some sort is involved. — In 
such cases no real significance is attached to the 
compressibility, which is averaged over the pressure 


interval containing the transition. In table 14 
corresponding data are given for the aluminous 


cements. 

On averaging the data for the dry hydrated 
portland cement specimens of tables 12 and 13 
regardless of age, water-cement ratio, or type of 
cement, an average compressibility of 3.1 107° 
atm”! is obtained below 5,000 atm and the corre- 
sponding value above 5,000 atm is 2.2 107° atm™'. 
The standard deviations of these figures are 0.7 


TABLE 13. Average compressibility of hydrated cement A 
containing different amounts of evaporable water 


water-cement ratio 0.5; l-week cure) 


Compressibility, 10-/atm 


Ey aporable 


water 1 to 5,000 5,000 to 10,000 
atm atm 
Percent 
32.4 14. 6 Ht 
26.9 12.8 5.7 
17.5 9.4 + 
5.2 3.1 2.5 
2 9 3.2 2.5 
0 2.5 2.2 


Average compressibility of hydrated aluminous 
cement of different age and evaporable-water content 


Water-cement ratio 0.5 
Compressibility, 10-° atm 
Curing period ~~ 
= 1 to 5.000 5,000 to 10,000 
itm atm 
Percent 

1 day . 0 15 1.8 
1 week 21.2 7.0 5 
Do 18. 1 5.8 4.1 

Do 15.8 18 3.5 

Do 0 3.0 29 

l month 0 7.0 3.7 
7's months 0 4.6 3.2 











10 atm and 0.3X10°° atm=', respectively. 
The single specimen of unhydrated cement powder 
had a compressibility of 2.810-° atm™' below 
5,000 atm and 1.6 107° atm above 5,000 atm. 

From the average compressibility and one other 


elastic constant, the elastic moduli may be cal- 
culated from the well-known relationships [5] 

k=3(1—2r)/B (9) 

(’=3(1—2v)/28(1+ 7), (10) 


where / is Young’s modulus, v Poisson's ratio, and 
G the shear modulus. Poisson’s ratios were deter- 
mined by sonic measurements on pastes made from 
cement A cured 1 week. They ranged from 0.27 
for a specimen of nominal water-cement ratio 0.25 
to a value of 0.31 for a specimen of nominal water- 
cement ratio 0.7. Adams and Williamson |28] have 
concluded that Poisson's ratio for plutonic rocks is 
approximately 0.27. These rocks are silicates and 
have essentially the same compressibility as hydrated 
portland cement. In table 15 elastic constants are 


TABLE 15 
from 


Elastic constants for the hydrated cements, estimate d 


compressibility measurements, Poisson's 


2; 


assuming 


ratio 


1 to 5,000 atm 5.000 to 10,000 atm 


$510) atm 6.310 atm 
Young's modulus 6.6 10° Ib/in. 9.310 Ib/in 2 
1.6105 kg/em 6.510 Kkg/em 
jt.S x10 itm 2.510 atm 
Modulus of shear 2.6 1 Ib/in. 3.7 10 Ib/in 
19X11 kg/em 2.610 kg/em 
[3.1 <1o-* atm 2.210 atm 
Compressibility 2.1107 in2/Ib 1510-5 in2/tb 
3.0X 10° em? /keg 2.1K10- em */keg 
o ’ . on . . 4 . 
given for dried cement pastes. The values for 


Young’s modulus and shear modulus have been 
computed from the average value of compressi- 
bility, assuming a Poisson’s ratio of 0.27. It may be 
noted that these elastic constants of cement are 
nearly the same order of magnitude as corresponding 
constants of the more compressible rocks, many of 
which are used as aggregates in concrete {5, 28]. 


The authors are indebted to E. T. Carlson at the 
Bureau for the specimen of tricalcium silicate, also 
to R. C. Valore, S. Spinner, and T. W. Reichart for 
the determinations of Poisson’s ratio by sonic 
methods, and to Vernon Dantzler for technical 
assistance. 


50 


(13] 


[14] 


H 


L.. 


S. 


Ir 
H 


C. 
G. 


F. 


C. 


R 
H 


* Birch, J. F. 


4. References 


_ L. Ritter and L. C. Drake, Ind. Eng. Chem. 17, 789 
(1945). 

_C. Drake and H. L. Ritter, Ind. Eng. Chem. 17, 787 
(1945). 

_C. Drake, Ind. Eng. Chem. 41, 780 (1949). 

* O. Slate, Comprehensive bibliography of cement and 
concrete 1925-1947, p. 160-162 (Joint Highway Re. 


search Project, Engineering Experiment Station, 
Purdue University, Lafayette, Ind.). 

Schairer, and H. C. Spicer, Handbook of 
physical constants, Geol. Soe. America Special Paper 
No. 36 (Jan. 31, 1942). 

_E. Weir, J. Research NBS 45, 468 (1950) RP2160, 
_W. Bridgman, Proc. Am. Acad. Arts Sei. 74, 11 (1940) 
. kK. Copeland and J. C. Hayes, ASTM Bul. No. 184, 
p. 70 (1953). 


Kk. Weir, J. Research NBS §2, 247 (1954) RP2496, 
_ E. Weir (unpublished data). 
. C. Powers and T. L. Brownyvard, Proc. Am. Concrete 
Inst. 43, 669-712 (1947): Bulletin No. 22, Portland 


Cement Assn. (March 1948). 

Kk. Copeland (personal communication from = manu- 
script Submitted to The Journal of the American 
Concrete Institute). 


T. C. Powers and T. L. Brownvard, Proce. Am. Concrete 


Inst. 43, 971-992 (1947). Bulletin No. 22, Portland 
Cement Association (March 1948). 

Giertz-Hedstrom, Proceedings of the symposium on 
the chemistry of cements (Stockholm 1938, published 
by the Royal Swedish Institute for Engineering Re- 


search , 


* C. Powers and T. L. Brownvard, Proce. Am. Concrete 


Inst. 43, 249 (1946): Bulletin No. 22, Portland Cement 

Assn. (March 1948). 

H. Adams, J. Am. Chem. Soc. 53, 3769 (1931 

Brunauer, P. H. Emmett, and E. Teller, J. Am 

Soc. 60, 309 (1938). 

iternational Critical Tables, IV 

Co., New York, N. Y., 1928). 

. F. Winterkorn, Soil Sei. 56, 109-115 (1943). 

i. Weir, J. Research NBS 49, 135 (1952) RP2349, 
Assarsson, Proceedings of the symposium on_ the 

chemistry of cements (Stockholm 1938, published by 

the Royal Swedish Institute for Engineering Re- 

search). 

M. Lea and C. H. Deseh, The chemistry of cement and 

econerete (Edward Arnold & Co., London, 1935). 

kK. Weir, J. Research NBS 64, 37 (1955) RP2562 

. H. Bogue, The chemistry of portland cement (Rein- 
hold Publishing Corp., New York, N. Y., 1947). 

. Lafuma, Rev. Gen. Sei. Appl. 1, 66 (1952), reference 

from Chem. Abstracts 46, 11617g. 


Chem 


(MeGraw-Hill Book 


. LeChatelier, Le Ciment 32, 82 (1927). 
Gessner, Kolloid Zeit. 46, 207 (1928): 47, 65, 160 
(1929). 
[28] L. H. Adams and E. D. Williamson, J. Franklin Inst, 


Wasuineton, May 3, 1955. 


195, 475 (1923 








7 


Jati 


f(r) 


equ 
wer 
pole 
pol: 
and 
mol 
inte 
of | 
effic 
frol 
‘| 
is ¢ 
Bes 
give 
res} 
the 
the 
suc 
era. 


Defe 








1, 782 


. 
4, 787 


t and 
\ Re. 


ation, 


rok of 
Paper 


60, 
940 
184, 


96, 


crete 





tland | 


anu 
rican 


crete 
tland 


ih On 
ished 


r Re. 


crete } 


ment 


hem 
Book 
19, 
the 
d by 
Re- 


and 





{ein- 


160 | 


Inst, 


Journal of Research of the National Bureau of Standards 


Vol. 56, No. 1 Research Paper 2647 


January 1956 


Formulas for Inverse Osculatory Interpolation 


Herbert E. 


Formulas for inverse osculatory interpolat 
formula. They cover the cases for n 


Salzer ' 


ion are obtained by inversion of Hermite’s 


2(1)7, where n is the number of points required in 


direct osculatory interpolation. The formulas provide an improved means for inverse inter- 
polation in the case where the first derivative is either tabulated alongside the function or 


is easily obtained. 


The author has previously given inverse interpo- | 
lation formulas [1, 2]? for finding r= sz 9+ ph from | 
f(z). in terms of f;=f(2,), where the 2#;=2 9+ 1h are | 
equally spaced at intervals of h. Those formulas 
were obtained by the inversion of Lagrange’s inter- 
polation formula. As the Hermite osculatory inter- 
polation formula [3, 4, 5] in terms of the function 
and its first derivative at preassigned points sz, is 
more accurate than the Lagrangian formula, inverse 
interpolation formulas obtained from the inversion 
of Hermite’s formula would be expected to be more 
efficient than the corresponding formulas obtained 
from Lagrange’s formula. 

There are many tables where the first derivative 
is either tabulated alongside the function (e. g., 
Bessel functions of the first or second kind, which | 
give J\(r) -Ji(r) |6, 7] or ¥4(x) Vo(r) [8], 
respectively, or probability functions [9] ), or where 
the derivative is very easy to obtain (e. g., tables of 
the more elementary functions and their integrals, 
such as tables of sine, cosine, or exponential inte- 
grals {10, 11] ). In all such tables, and in many 








others, the user will find the inverse interpolation 
formulas given below to be particularly convenient, 
especially when the tabular interval / is too large for 
sufficiently accurate inverse interpolation, using the 
formulas in terms of the functional values alone. 

These present formulas give p in terms of f(s9- 
pry=f, fled=f,, and f’(a)=fi for 25=29+th at 
equally spaced intervals 4, where the 7 ranges from 
—| (n-1)/2] to [n/2], for n=2(1)7, where n is the 
number of points required in direct osculatory 
interpolation. Although the direct interpolation 
formula for n=6 and n=7 is of the 11th- and 13th- 
degree accuracy, respectively, the inversion formula 
for p that is given below does not go beyond the 
10th- degree terms (and in most practical problems 
it is verv rarely that one will go that far; in fact, 
the first few terms will usually suffice). Thus we 
do not use the coefficient of p'' in the 6-point Hermite 
formula nor the coefficients of p"', p'?, and p* in 
the 7-point Hermite formula. For every n, we 
define r= (f—fo)/hfg, and corresponding to n=2(1)7, 
quantities s, t, u,v, w, 4, y, 2, and @ are defined 
as follows: 


a=2 
s=(—3{ fo—fi} —h{2fot+Si} hho, 
t=(2{ fo—fhi} +hi fot Si} /hto, 
u=r=w=r=y=7=7=0. 
n= 
e=(fi—2fo fee Se ii\) hfo, 
4 | | 


4 | 4] 
1 | . 
u (—5 ‘f l 2hort hi 
3S | : h\. 
/ G I l hry iJ 
w= 2=y= 2=2=0. 


Present address, Diamond Ordnance Fuze Laboratory, Department of 
Defense, Washington, D.C 
- Figures in brackets indicate the literature reference at the end of this paper 
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if T Sf; T ty ) hf, 


Fes fi, ) hts, 


f'..+4fe4 fi ) hfy, 








t=(— 


n 
| ) 
( 56 Ff 9Q7F 216¢ one! h »f’ 4 + r| ol 
108 )o »f j~ —* ‘fi | 216/; T 25f> T79 2) 1 ISS, LSJ, - fa) ) hfo, 
t ( mes 124/ 7 f, —108f,+11fo;— A Sf 634; i" hf; 
108 |} °~* _ ' /2! 36 ee a8 ) So, 
] = ° 2 » h = -s - | “9 
u las 25(F_y / 139(fo—h) +36 f_,+72f)+ 63), 4fz, ) hfo, 
( 59f_,—54f.—27h, +22 fo; + . f" 1sf.—9f,—2 | ) hf; 
108 _— ogilionss ” 36 |*° =! —— "JO 
l , h 
; oe —_ =“) ‘ ) —— - cr’ J <« — ‘ id 7 
r=( aa 152 $1 fy—108f,—25 fs} ge | BS + 36fo+ 27K, + 2fy ) hf. 
l — h 
. | s “ ; 
A pe iC fy) +27 (fof) tag Sat Pot Mit fe ) hfi, 
Y = () 
n=5 
a ap i h ‘ 
- a ( f_2+hfs 128(f / 2T70f, +35 (f f3)+32(f-,—S;) ) hf; 
Cn l a1({ f.) 4+ 1408 . h 2 , ee vieaiial “ 
864 |° J2 J Ji) “144 (J + fo) + 64(f f,)+ 3607, ) hf 
] — ; . ; h , — - ss 
( - “FRR 41( fos 256(f f,)—594f, - 30 f,)+64(f_,—f) ) Afi. 
1152 a— J 92( 7 f,) T 99 3(f f,)+128(f_,+f)) d9OT 4 ) hf. 
I — , ' - ; 
t > 13(f-2+ fe) +32(f_1+f,) —90fo; + 3 (fF) +1208. —S1) ) hf. 
(— 29(f f 7 2(f an fs 4 (f OF, p 
576 oo 06 | lt 11) + OOJo ) hf, 
(— : ll( fio+f L6( 7 547 h F. fs _ 7 c 
Sti4 a ‘ h)— 94 fo Seq | S-2—J2) T8S-1— Si) ) hf, 
] ” ; . h 
» . P er .? . ad -* ava Ue | “7 
Fa 25(f2—fx)+160(f—f)) tare | Sat S24 160+ Si) 36/0) ) hf 
(0. 
n—6 
aa . 
DOD Vf 37 Of S500! SOOO! LOO0/, + 491 
B® low sean ee - 
ton 3 150f"., —400f, —600f; — 75/3 —2f; ) hf, 
I esewe eee 
108000 | oem 2 LO9S00S SO00/, — 168000/, — 2625/, + SASS, 
Seng ) lf -2t 1900/ + 8600/0 + 1200f, — 75f ifs) ) hf. 
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l ; — iit aat ae ‘ : phere is , 
u ( 51600 943f_.— 2050f_ , —62550f, + 45200F, + 175751, + 882, 
h re aa ae on . ‘ 
7440 | 132+ 260/~, — 2400/0 -- 2960; — 445f;— 12/3 )/ afi, 
(nes00 5729f_.+121400f_,— 16000f,— 118400f, 4 6175s, + 1096f, 
ao 
, h | oo af = mT) o 7 7 , , 7 
+o (3éf 4+ 1720f_, + 5140f, + 1120f, — 155f, 8/5) ) hf, 
| , ; ; 
we ( 144000 ) 1274f_.--S84875f_,— 189000/,+ 175000/, + 92750f,4 4851f, 
h / on 7) 7 af , ot os 
+ 74490 | Dbl 2 2175!" — 19800/, — 20100f; — 3600/; 99/5) ) hf, 
| = . _ , - dca 
r (“7a 5743f_.4 50125f_,— 22000f, — 46000f, = 10625/,+ 1507f; 
, h 9 9 P | 
———_ 1 30f",- 925f 1900F) + 100; — 225f; 11fs) ) hf; 
4800 | v2 ae 
| ‘iit 9 tiles mbieee eo — 
Yu licen 172f_.+3275f_,+3900/,--4600/, — 2600/,— 147f; 
ee , 
a f’.+45¢'.,4-200f.+ 180° +35f2+f2! ) r? 
| $80 2 4 ) 1 f ] ry 5 fo Js, hf, 
: ( | 1 351f_.4+-750f ,—4000f)-+2625/,+ 274 
: 86400 _— — Pe eave ee ae 
h asin a - ia ~ - 
+5959 ) Bf 2+ 8Of@ 1 — 20f— 120f; — 45), 2f3) ) hf, 
= | ,; = : ge ee . | 
: (— 5 '1666f_.+ 10625f_, + 9000f,—13000/, —7250/.— 441 fy. 
$2000 . Js 
t -/ =~ pe > 7 ~ 77 7 ort -r 
-—vT Sf. +175f,+ 600f) + 500f; + 100f5 4+ 3/; ) hfi. 


n=7Z 


l ; : ; ) 
57 (f_3+fs) +4968 (f_2+-f2) +-43875(f_14 9S000f,, 
( saaae 157 (f_3+fs) +4968 f_2+-fo) +-438875(f_1 +f) 000f 
i) ; | 


fs fa) +54 Sf) + 075 Ff), ) Ai, 


T7200 | § 





le fe _ 4) 
(— {98000 167(f_3—fs) + 8667 f_2—f2) + 192375 (f_1—fi) 
ih = 


-s - - P } 
_ } 1. F. bat o-+-f.)-+- 2005 Ley 8 op | - 
3600 | f-3 tf) + 810-2 + fe) + 2025 F_ 1+ fi) 9800/,, ) hf. 
l : : ) 
— 73: 19] 2726 ! | ga11988 pei cae 
u ( 648000 7339(f_3+f3) + 213786 (f_24+-fe) + 631125 (f_, +f) — 1704500f, 
ae ) 
51600 |47 F-2 fa) +2403 F~ 2—f2) + 20925(F fi);) hfi,, 
l Ps . "a ° , - nia mia ‘ 
( 1296000 5206 a—I3)-+4 24983 1 f 2 fo) L 9517 750 ' fy ) 
h “g ; ) 


94(f@ s+ fa) + 7209(F" 2 +f) + 125550(f~ +1) +340900f5) )/ hifi, 
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: 129600 | 

| I ms) or ° . 7 
} 5P109(f_ a4 )+- 1B 10016( +f.) + 15258751 | varie | 

u ( 5184000 ) fests (f_oth, 29879 (f_14-4,)—57 76000fo 


h Jagr ¢ , - 7 7 a wel 
+ 775800 }3370 s— fi) + 15648 (f~2—f3) + 88275 (1 — fi); ) hfs, 
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18493 (/ 


z==( fy 


772368 (f 








ou ) ) 
5184000 | 
h - - —— i ' Pe , 1 
—- 337(f fo 23472(7 +-f5)+264825(f_1+7,)+577600F, ) hf; 
518400 : : | | 7 | 
] ; j = ; ipa al } 
y (— ‘. 2174 f 117408(1 f.)+72375(f f 392000F 
. 1728000 ; ’ ; 
h » i . ; e sc acl ‘ : 
omen f’ s—f3) + 1584 (ff) + 6675-1. —f1)} ) fe, 
5/6000 ' : sted ' . 
1 ; =(\0 ; 979GOS 
(— 2209 (f_;—f,) + 70584 (f_»—fe) +272625(f_:—f; 
1728000 ; , } 
h , , icc i . a - ? 
a f f°) +2376(f'.o+f 20025(f., +f) +39200f, ) Af’. 
172800 : =" , : “hg 
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